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Preface 
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Abstract 

This  thesis  is  an  evaluation  of  the  factors  that  determine  the 
maximum  acquisition  and  tracking  (lock-on)  range  of  laser-guided 
systems  in  the  military  close-air-support  mission.  The  problem  is 
divided  into  a  study  of  the  designed  system  parameters  which  are 
characterized  by  a  clear-air  lock-on  range,  and  a  study  of  the  effect 
of  the  intervening  atmosphere.  Three  model  atmospheres  are  considered: 
a  homogeneous,  an  exponential,  and  a  three-layer  stratified  atmospheric 
aerosol  vertical  structure.  The  vertical  aerosol  number  density 
distribution  of  a  real  atmosphere  can  be  bracketed  between  a  homogeneous 
and  an  exponential  atmospheric  model.  Conversion  charts  from  clear-air 
lock-on  range  to  actual  lock-on  range  for  the  three  model  atmospheres 
are  plotted.  It  was  found  that  beamspreading  by  aerosols  can  be  large 
for  moderate  attenuation  coefficients.  Reflective  cross  sections  of  a 
jeep  model  were  found  to  have  a  standard  deviation  of  +20%  of  the  mean 
for  two  olive  dr~b  paints.  Attenuation  coefficients  at  1.065  microns 
are  determined  for  atmospheric  gases  and  rain.  A  method  is  presented 
for  determining  aerosol  attenuation  at  1.065  microns  from  the 
visibilities  measured  through  red  and  blue  filters;  comparison  of  this 
method  with  the  computer  analysis  by  others  was  successful  but 
experimental  tests  were  inconclusive  due  to  equipment  and  weather 
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LOCK-ON  RANGES  OF  LASER-GUIDED  SYSTEMS 
I.  Introduction 

Background 

The  purpose  of  this  study  is  to  investigate  the  factors  that 
determine  the  lock-on  range  of  laser-guided  systems  in  the  military 
close-air-support  mission. 


A  laser  beam  is  used  to  illuminate  a  target  as  in  Figure  1.  The 
reflected  laser  energy  is  detected  by  a  narrowband  electro-optical 
receiver  and  used  to  guide  the  receiver  to  the  target.  The  maximum 
range  at  which  the  receiver  can  acquire  end  track  the  target  (lock-on 
range)  is  limited  by  the  laser  power,  beam  divergence,  target  reflective 
characteristics,  receiver  sensitivity,  and  the  atmosphere  in  the  path. 

Atmospheric  attenuation  is  a  severely  limiting  factor  in  the 
operational  employment  of  laser-guided  systems,  since  the  usefu]  lock-on 
range  may  bo  reduced  by  a  factor  of  20  in  weather  that  a  pilot  car.  still 
fly  by  visual  references.  Certain  laser  wavelengths  such  as  1.065 
microns  (4)  have  better  propagation  characteristics  than  those  in  the 
visual  range,  due  to  less  scattering  and  absorption.  In  adverse  weather, 


1 


GEP/PH/72-12 


attenuation  at  1.065  (i  is  principally  due  to  scattering  by  liquid  water 
aerosols,  plus  slight  absorption  from  water  vapor  and  from  dust  aerosols. 
Prediction  of  attenuation  in  a  real  atmosphere  is  beyond  the  state  of 
the  art,  as  is  prediction  of  the  reflective  characteristics  of  a 
complex  target  in  a  combat  environment.  Finally,  there  is  very  little 
standardization  in  the  literature  regarding  laser  beam  divergence 
definitions,  reflectance  terminology  and  characteristics,  and 
atmospheric  transmittance  terminology  and  characteristics. 

The  solution  of  this  problem  is  of  interest  to  the  Operations 
Evaluation  Group  (USAF-SAV)  of  the  Air  Staff.  In  the  presentation  of 
this  problem,  they  instructed  the  researcher  to: 

1.  Determine  what  factors  most  severely  limit  the  maximum  lock-on  range 
of  laser-guided  systems  in  adverse  weather. 

2.  Propose  improvements  to  increase  the  lock-on  range. 

3.  Comment  on  factors  that  may  affect  the  tactics  or  operational 
employment  of  these  systems,  but  keep  the  written  thesis  unclassified. 
Scope 

Important  sub-problems  are  beam  characteristics,  beam-target 
geometry,  target  reflectance,  receiver  characteristics,  attenuating 
mechanisms,  aerosol  compos i  ion  and  size  distribution,  atmospheric 
vertical  profile  models,  real  atmospheres,  statistical  variation  in 
lock-on  ranges,  and  methods  of  approx 'mating  values  of  aerosol 
attenuation.  The  tasks  are  to  determine  the  state  of  the  art,  identify 
areas  where  knowledge  is  lacking,  evaluate  the  relative  importance  of 
these  areas,  propose  solutions,  and  test  these  proposals  experimentally 
or  by  computet  analysis. 
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Limi tat: ions  and  Assumptions 

This  thesis  is  limited  to  an  analysis  cf  1 . 065— (JL  laser  systems  at 
ranges  of  0.5  -  10  Km.  Nonlinear  effects  and  optical  countermeasures 
are  not  considered.  It  is  assumed  that  beam  wander  is  either 
compensated  for  in  laser  pointing  and  tracking  or  is  negligible.  It  is 
assumed  that  the  human  (photopic)  eye  can  detect  a  2/  contrast  through 
narrowband  filters  anywhere  in  the  visual  range  as  well  as  with  white 
light. 

Organization 

The  clear-air  lock-on  range  due  to  system  geometry  will  be 
determined  in  Chapter  Two.  Beam,  target,  and  receiver  characteristics 
will  be  discussed  as  they  apply  to  the  operational  problem.  Chapter 
Three  will  evaluate  the  attenuation  of  laser  beams  by  gases  and 
aerosols  in  model  and  real  atmospheres.  A  method  for  determining 
attenuation  coefficients  from  meteorological  measurements  will  be 
presented  and  evaluated.  In  Chapter  Four  the  clear-air  lock-on  range 
will  be  combined  with  three  different  model  atmospheres  to  yield  a  true 
lock-on  range.  The  statistical  variation  in  this  true  lock-on  range 
will  be  discussed.  Chapter  Five  will  test  and  evaluate  one  of  the 
proposals  to  determine  the  attenuation  coefficients  from  meteorological 
range  measurements  through  narrowband  filters.  The  major  conclusions 
and  recommendations  are  given  in  Chapter  Six. 
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II.  Geometry  of  the  Illuminator-Target-Receiver  Problem 
Characteristics  of  the  Laser  Beam 

Military  lasei  specifications  require  a  beam  with  a  gaussian 
intensity  profile  .--.cross  the  beam.  This  may  be  achieved  in  systems  with 
electrically  polarizing  cells  as  the  Q-switch,  but  is  almost  impossible 
to  achieve  in" systems  that  use  a  rotating  mirror  as  the  Q-switch.  The 
latter  device  almost  always  results  in  a  multi-mode  lasing  action,  with 
highly  irregular  and  erratic  intensity  distribution  across  the  beam.  The 
difference  in  beam  pattern  is  mostly  of  academic  and  computational 
interest,  as  scintillation  in  the  atmosphere  will  often  cause  far  more 
serious  disturbances  to  the  beam’s  internal  structure. 


The  computational  problems  involved  in  a  multi-mode  beam  can  be 
immediately  understood  by  observing  Figure  2,  a  series  of  cross¬ 
sect  'onal  sweeps  across  a  military  1 . 065— {j.  laser  beam  by  the  Target 
Signatures  Analysis  Center  (TCAC),  University  of  Michigan  (Ref  2:51). 

If  several  pulses  re  averaged  together,  the  irregularities  in  the 
intensity  dii,  ibution  average  out,  so  for  many  practical  cases  an 
assumption  of  .onstant  illumination  across  the  beam  is  better  than  an 
assumption  of  a  gaussian  distribution.  This  report  will  consider  two 
cases  of  laser  beam  profiles:  a  gaussian  and  a  constant  average 
intensity  (irradiance)  distribution  across  the  beam,  thus  covering 
most  practical  cases.  Detection  warning  and  countermeasures  receivers, 
however,  can  respond  to  the  peak  pulse  intensity. 

The  constant  average  illumination  beam  is  defined  by: 

I(r)  =  I(o)  for  r  <  w 

CD 

I(r)  =  0  for  r  >  w 

The  general  definition  of  a  gaussian  intensity  distribution  is: 
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Kr)  _ 
I(o) 


(2) 


where: 

r  =  radius  from  the  cenf.er  of  the  beam  pattern,  in  cm. 

w  =  a  defined  hali— width  of  the  beam  pattern,  in  cm. 

c  =  a  scale  constant  depending  on  the  definition  of  w  (see  Table  I). 

c  =:  2.3026  for  the  standard  definition  used  in  this  thesis. 

I(r)  =  irradiance  at  radius  r  in  the  beam  pattern,  in  watts/cm 

I(o)  =  irradiance  at  the  center  of  the  beam  pattern,  in  watts/cm“. 

c  f  I(t ) 

— ^  =  a  height-normalizing  constant  such  that  j  ~  dS  =  1. 

ttw“  *3 

S  “  integrating  surface  across  the  beam  through  r,w. 

There  are  several  means  of  specifying  whole-angle  beam  divergences 

of  a  gaussian  laser  beam  in  the  optical  far  field.  All  appear  with 

nearly  equal  frequency  in  military  laser  illuminator  specifications 

without  much  consistency  or  clarification  about  which  method  is  used. 

In  more  recent  reports  and  specifications  the  beam  divergence  is  most 

commonly  defined  as  the  cone  containing  9 0%  of  the  total  beam  energy. 

This  thesis  will  follow  this  convention.  Of  all  the  definitions  used, 

this  results  in  the  narrowest  real  beam  for  a  given  numerical  divergence. 

Beam  divergences  are  sometimes  defined  by  the  cone  containing  50% 

of  the  beam  energy;  or,  by  the  cone  at  the  points  where  the  local 

2 

intensity  becomes  10%,  50%,  1/e,  1/e  ,  -10  dB,  or  -3  dB  of  the  intensity 
at  the  beam  axis.  Caution  must  bo  used  to  determine  which  definition  of 
beam  divergence  is  used  with  which  laser  illuminator.  A  computer 
program  was  developed  (Appendix  A)  to  compare  all  other  definitions  with 
the  first.  The  results  are  in  Table  J. 

Irradiance  of  the  Target 

Irradiance  by  the  Laser  Beam.  Once  the  laser  beam  has  reached  the 
optical  far  field,  the  beam  diameter  can  then  bo  determined  in  the 
absence  of  atmospheric  effects  b, . 
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d  =  2  R  Tan  (^)  +  d 
l  o 


(3) 


where: 

a  =  beam  divergence  whole-angle  in  Che  far  field# 

R  =  illuminator  -  target  range, 
d  =  beam  diameter  m  the  far  field  at  range  R. 

d  =  laser  illuminator  beam  exit  diameter  (if  this  is  in  the  far  field) 

°  or  the  diameter  of  the  rearward  projection  of  the  divergence  angle 
to  the  laser  exit  window. 

If  a  <  100  milliradians,  Eq  (3)  is  approximated  within  0.4%  by: 

d  =  Ra  +  d  (d,d  in  cm;  R  in  meters;  a  in  milliradians)  (4) 

To  0  0 

The  power  density  (irradiance)  in  the  beam  is  of  considerable  value 
in  optical  countermeasures  design.  Two  cases  are  of  primary  interest: 
the  constant  average  distribution  and  the  gaussian  distribution.  In  the 
absence  of  atmospheric  effects,  the  former  is  given  by: 


P  =  transmitter  power. 

d  =  beam  diameter  enclosing  100%  of  the  constant  beam  energy. 

H  =  irradiance  at  the  target  position  in  the  beam, 
a  =  whole-angle  beam  divergence  at  100%  enclosed-energy  points. 
dQ=  beam  exit  diameter,  as  in  Eq  (3). 

For  a  gaussian  distribution,  the  irradiance  at  the  target  position  is: 


(from  Eq  (2)): 


hcd  oisscger  e-2-302<t-fy] 

TT  (?)  V 


-  9. 2104 P e~9. 2104(4) 

2  a 
TTd 


where : 


d  =  beam  diameter  enclosing  90%  of  gaussian  beam  emergy  =  Ra  +  d  . 
a  =  beam  divergence  whole-angle  between  90%  enclosed-energy  points, 
r  =  radial  position  normal  to  the  beam  at  range  R. 

P,  dQ,  H,  R  as  in  Eq  (5). 

Irradiance  bv  Other  Sources.  Of  all  the  other  sources  of 
illumination  of  the  target,  the  direct  solar  radiation  is  by  far  the 
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most  intense.  For  a  100  k  window  about  1.065  p,  at  sea  level  on  a  clear 

2 

day  it  is  3  -  6  watts/meter  ,  depending  on  the  solar  zenith  angle  (Ref 

37:  6-6, 6-8).  The'  average  radiance  of  the  scattered  light  from  each 

**6 

corresponding  solid  angle  of  the  clear  daytime  sky  is  2x10  as  strong 

as  the  radiance  from  the  sun,  at  1.065  p  (Ref  3:1314),  but  over  the 

whole  sky  the  total  becomes  21%  as  much  as  the  direct  solar  irradiance. 

Individual  bright  clouds  may  have  a  radiance  up  to  ten  times  that  of  the 

background  sky,  but  thicker,  darker  clouds  may  greatly  attenuate  the 

direct  and  indirect  irradiance  from  the  sun,  so  that  the  net  contribution 

from  clouds  is  unpredictable:  it  varies  with  the  sun  angle  from  the 

target  and  the  angle  from  the  target  to  each  of  the  individual  clouds. 

The  total  day  or  night  emission  from  sky  or  clouds  at  1.065  p,  is 
-3 

less  than  10  as  strong  as  the  sky  or  cloud  scattered  light  (Ref  3: 

1314),  and  can  be  neglected.  Thus  the  total  daytime  direct  and  indirect 

2 

background  irradiance  onto  the  target  is  3  -  8  watts/meter  (300  -  800 
2  o 

microwatts/cm  )  for  a  100  A  window  at  1.065  p.  This  is  a  factor  of  10^ 

g 

to  10  lower  than  the  peak  pulse  illumination  by  most  typical  systems 
at  reasonable  target  ranges. 

Portion  of  the  Beam  Energy  Intercepted  bv  the  Target 

Since  the  beam  at  the  target  may  be  smaller  or  larger  than  the 
projected  cross  section  of  the  target,  it  is  necessary  to  compute  what 
portion  of  the  total  beam  energy  is  intercepted  by  the  target.  The 
remainder  will  be  lost  into  the  background,  which  usually  has  a  higher 
reflectance  and  may  bo  a  significant  distance  behind  the  target.  This 
results  in  two  separated  illumination  and  lock~on  problems. 

A  computer  program  was  created  to  determine  the  portion  of  the 
laser  beam  intercepted  by  targe ‘■'s  of  different  shapes:  round,  square. 
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.and  rectangular  for  length-to-haight  ratios  of  n  =  1.5,  2,  3,  4,  5,  6, 

8,  10,  and  20j  for  a  gaussian  beam  and  a  constant  cross-section  beam. 

The  irradiance  of  Lqs  (5)  and  (6)  was  integrated  across  the  target  shape 
and  the  results  pljtted  in  Figure  3  and  Figure  4.  The  calculations  are 
in  Appendix  B.  In  figure  4  there  was  no  visible  change  in  the  curve 
shape  above  n  =  8. 

Target  Reflectance 

A  real  surface  can  be  considered  to  be  covered  with  a  broad  size 
distribution  of  Mie  scattering  particles.  The’ surface  will  have  a 
relatively  broad  and  strong  forward  scattering  lobe  (specular  beam)  and 
a  general  diffuse  reflection  in  all  directions,  a  result  of  summing  the 
Mie  reflection  lobes  of  all  the  particles.  The  real  surface  will  also 
have  interference  effects  due  to  the  surface  fine  structure  (grain  of 
wood,  cloth)  and  polarization  effects. 

The  term  "reflectance"  is  frequently  used  synonymously  with  albedo, 
spectral  reflectance,  directional  reflectance,  bi-directional 
reflectance,  reflection  distribution  function,  and  target  reflectance 
cross  section.  It  is  so  easy  to  misuse  the  published  numerical  values 
of  these  terms  that  a  few  words  of  caution  and  clarification  are  in  order. 

Albedo  is  usually  the  same  as  reflectance,  the  ratio  of  total 


radiant  power  reflected  P^  to  the  total  power  incident  P^  on  a  surface: 


But  albedo  is  sometimes  used  to  include  black-body  emission.  Spectral 
reflectance  is  the  reflectance  in  a  narrow  band  about  a  wavelength  X: 


(8) 
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Directional  Reflectance  is  the  ratio  of  the  radiant  power  reflected  into 
a  specific  solid  angle,  to  the  total  power  incident  from  all  directions: 


P_(e  . 0  )  j(e  ,0  )  dQ 

pd  =  P<W  z  =  — f-E - 

i  i 


for  a  body 


(9) 


where: 


**r^r,0r^  power  reflected  into  direction  (6^,0^),  in  watts. 

J(V0r}  *  radiant  intensity  into  direction  (6^),  in  watts/steradian. 
dQ  =  solid  angle  subtended  by  the  small,  distant  receiver. 

For  a  small  target  surface  area  dA  ,  this  becomes: 

Nr(V0r)  dQdA^  dQ  cos  e 


p,  =  p(e  ,0  )  =  — — E— r 

Md  -  MV  r>VTJ  H  dA 


H 


Cio) 


where : 

H  =  irradiance  onto  the  actual  target  area  dA  in  watts/cm2  of  actual 
target  area.  ,  p  t  ~ 

Nr^r,0r^  -  - 2 — -p — ~  radiance  in  watts /steradian-cm2  (n) 

cos  dAt  df)  of  projected  area  perpendicular  (Aj. ) 

to  the  line  to  the  receiver. 

A  very  dusty  surface  is  approximated  reasonably  uoll  by  a  perfectly 

diffusing  (Lambertian)  surface  =  constant  =  si),  so  the  total 

power  reflected  from  a  target  surface  area  dA  is: 

t 


dP 

dA 


E 

2 TT  2 


>  t-u  i 

T  =/n  cos  ^  dQ=/  /n  cos  6  sin  6  d6  d0  =  Nit 

t  Q  i  J0 

For  this  very  dusty  surface  the  directional  reflectance  is: 


(12) 


dP 

».-X 


dP. 

_ i 

dA_ 


-  Eii  xt  H 

-  *tt—  ,  or  N  - - 

H  tt 


.  9 

m  watts/steradian**cm“  (13) 
of  projected  area  perpendicular 
to  the  line  to  the  receiver. 


This  "projected  area"  is  either  the  projection  of  the  illuminated  target 
surface  area  onto  a  surface  perpendicular  to  the  line  to  the  receiver 
CAt  cos  er),  or  the  projected  area  that  is  within  the  field  of  view  of 
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the  receiver  at  that  projected  point,  whichever  is  least.  The  irradiance 

H  is  the  power  onto  the  surface  divided  by  the  actual  target  surface 

2 

area  A„,  in  watts/cm  .  It  must  be  uniform  and  cover  all  of  A.  ,  or  A 
t  t  c. 

must  be  broken  up  into  smaller  elements.  For  a  real,  uon-Lanbertian 
surface  N  usually  has  a  very  strong  angular  dependence  other  than  the 
cos  6^  dependence,  and  the  directional  reflectance  of  the  surface 


*  if/«V 


0  )  cos  8  sin  8  d6  d0  becomes  very  difficult  f  n 


r  r  r 


evaluate.  pd  may  be  much  greater  than  1  at  certain  angles  but  will 
average  out  over  all  space  to  be  1  or  loss. 

Bi-directioiial  reflectance  (partial  reflectance)  (reflectance 


distribution  function)  is  the  ratio  of  the  radiant  power  reflected  into 

a  specific  solid  angle,  to  the  radiant  power  incident  from  a  specific 

solid  angle.  It  is  assumed  that  the  incident  power  is  collimated.  The 

surface  irradiance  is  therefore  H(8.,0.)  =  H,  cos  6..  Similarly  to 

the  arguments  for  directional  reflectance,  there  is  for  this  geometry 

a  bi-directional  reflectance  p  =  p(6->0-j8  ,0.)  such  that  C  d Q=  p  , 

p  lirx  */'P  d 

where  fi  is  the  influx  or  efflux  hemisphere  above  the  target  surface. 

For  a  perfect  Lambertian  surface  this  becomes  (Ref  27:117): 


pi  =  pP/dn 


Therefore  N  =  becomes  N  =  p  H.  For  a  real,  non-Lambertian  surface: 
if  P 


J  N(8  ,0  )  cos  8  sin  8  d6  d0 
jj  '  r  r  r  r  r  r 


H(8.,0.) 


where  the  integration  limits  are  over  the  reflectance  cone  d0^  d0^. 

For  a  real  surface  p^  may  bn  much  greater  than  ~  at  certain  angles 
but  will  average  over  all  space  to  be  ~  or  less.  Integ  ation  of  the 
receiver  or  source  flux  over  the  total  hemisphere  above  the  surface 
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yields  the  total  reflected  or  incident  radiation,  respectively,  and 
hence  yields  the  directional  reflectance  which  is  tt  times  the  average 
ppi  Pjj  =  IT  pp  •  Integration  of  both  source  and  receiver  fluxes  over  the 
hemisphere  above  the  surface  yields  the  total  reflectance  p: 


P  =  Pd  58  *  Pp  <16) 

The  term  "directional  reflectance"  is  frequently  and  unfortunately  used 
when  an  author  means  "bi-directional  reflectance".  The  "bi"  is  implied. 
Reflective  Cross  Section 

There  are  two  parts  to  the  problem:  target  material  and  target 
geometry.  The  Target  Signatures  Analysis  Center  has  measured,  computed, 
or  catalogued  the  spectral,  directional,  and  bi-directional  reflectances 
of  several  thousand  materials  of  tactical  interest.  Their  mathematical 
model  for  the  bi-directional  reflectance  of  a  material  consists  of  a 
weighted  sum  of  three  parts:  one  part  perfectly  diffuse  reflection  from 
a  hard  surface,  one  part  perfectly  specular  reflection  from  a  hard 
surface,  and  one  part  multiple  diffuse  volume  reflection  from  a  thin 
volume  of  scattering  dielectrics  on  the  surface  (paint).  The  three 
parts  are  weighted  according  to  theoretical  and  experimental  data. 

For  a  complex-shaped  target  such  as  a  vehicle,  the  different 
surface  materials  at  different  angles  results  in  a  complicated  total- 
target  bi-directional  cross  section.  For  constant  irradiance  H: 


°c5i4/pP,!t,s 


m 


watts /stored ran 
watts /cm2 


(17) 


^Total  vehicle  surface 


TSAC  has  computed  this  function  for  a  mathematical  jeep  model 
which  was  validated  by  field  measurements.  Oik*  set  of  their  computations 
is  in  Figure  5,  a  polar  logarithmic  plot  of  the  total-target  reflective 


cross  section  a  for  three  receiver  xenith  angles.  The  fj.eculai  peak 

U 
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PER  IRRAO  IANCE  OP  1  N  A  T  T  /  C  ' 

Fie.'  6.  Statistical  Variation  in  Reflective  Ctoa.'  Sect:.*,  : 
Random  Annular  Orientation 


probability  of  reflective  cross  section  being  above  each  range 
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PER  IRRflGIANCE  OF  1  WflTT/CM* 

Fig.  7.  Statistical  Prob.dii  1 :  ty  of  Jeep  Having  Reflective  Cross  Section 
Above  the  Indicated  Value 
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from  the  flat  vehicle  top  is  conspicuous.  Samples  taken  from  the 
jeep  used  in  the  validation  tests  had  reflectances  at  1.065  |i  for  new 
paint  of  2.9- 3.3%,  old  paint  of  4. 8-6. 7%,  black  canvas  of  9.3%,  and 
olive  drab  canvas  of  12.1%  (Rtf  35).  Most  natural  vegetation  and  dirts 
have  a  higher  reflectance  range:  10-60%. 

Computed  values  of  c^.  were  measured  from  TSAC  charts  (Ref  2,  Ref  34) 
and  the  statistical  probability  of  getting  a  certain  reflective  cross 
section  for  completely  random  angular  orientation  was  computed  and 
plotted  in  Figure  6.  The  receiver  was  unpolarized;  the  transmitter  was 
either  horizontally  or  vertically  polarized.  The  probability  of  the  o 
being  above  any  value  was  computed  and  plotted  in  Figure  7.  The  bottom 
graph  of  Figures  6  and  7  is  of  a  very  specular  paint,  not  likely  to  be 
encountered  in  tactical  combat  environments,  but  is  presented  for 
comparison  to  the  military  olive  drab,  low  reflectance  paint,  on  the 
same  jeep  model.  These  computations  are  for  generally  clean  surfaces 
(not  mud-caked,  oily,  or  dusty).  A  dusty  surface  has  generally  a  much 
more  diffuse  character  and  a  much  higher  reflectance  (10-60%).  Wot 
surfaces  are  completely  unpredictable.  This  procedure  for  the  analysis 
of  reflective  statistics  of  potential  threat  vehicles  yields  the 
maximum,  minimum,  and  distribution  of  for  targets  under  clean 
conditions.  It  is  strongly  recommended  that  these  statistics  be 
computed  for  many  enemy  vehicles  of  potential  threat  value,  under  many 
different  field  conditions,  for  wet  and/or  dirty  surfaces. 

Lock-on  Range 

Of  the  energy  lc'lccted  by  the  target  in  a  general  direction 
towards  the  receiver,  J(®r»^r)  watts /stcradian,  only  the  portion 
entering  the  effective  receiving  area  will  be  available  for 
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amplification,  detection,  and  electronic  analysis.  The  detector  of 
effective  area  (real  area  decreased  for  losses  in  receiver  optics) 
and  detector  -  target  range  subtends  a  solid  angle  from  the  target  of: 

df)  =  d  - d  (is) 

Ed 

where  6^  is  the  angle  between  the  detector  normal  and  the  line  to 
the  apparent  target  through  the  receiver  optics.  The  cos  is  necessary 
because  the  detector  may  not  be  pointed  along  the  line  to  the  target 
before  lock-on. 

If  the  illuminating  beam  is  larger  than  the  target  and  has  uniform 
irradiance,  then  for  zero  atmospheric  attenuation  the  power  received  by 
the  detector  in  watts  becomes  (from  Eqs  (5),  (17),  (13)): 

A PA  .  a.  cos  6 , 

P  =  Jd()  =  H  a  dQ  =  — 5 - y  in  watts  (19) 

ttrI  (Ra  +  d  ) 

d  0 


If  Pdm  is  the  minimum  detector  power  in  watts  necessary  for  the  receiver 
to  acquire  and  track  (lock-on)  the  target,  then  the  maximum  lock-on 


range  is  given  by  (from  Eq  (19)  (for  zero  atmospheric  attenuation)): 


(20) 


If  the  irradiance  onto  the  target  is  not  constant,  then  the  energy 

distribution  inside  the  beam  has  to  be  included  in  a  new  reflective 

*  -1 

cross  section  c  (with  units  of  steradian  ): 

°t  =  P  =  p  /  PpI!<VV  dS  (steradian-1)  (21) 

S 


If  the  illuminating  beam  is  not  much  larger  than  the  target  and  can  be 
approximated  by  a  gaussian  or  constant  irradiance  beam,  then  Figure  3  or 
A  will  yield  the  portion  K  of  the*  total  beam  energy  intercepted  by  the 

t 

target.  Then  J  -  PKcJ^.  and  Eq  (19)  becomes: 
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The  maximum,  minimum,  and  distribution  of  lock-on  range  can  be 
determined  by  substituting  the  appropriate  values  of  0^  or  c  from 
Figures  6  and  7.  This  concept  of  a  clear-air  lock-on  range,  independent 
of  atmospheric  conditions,  and  fully  descriptive  of  system  geometry  and 
its  statistical  variations,  will  be  used  frequently  in  this  thesis. 

It  is  generally  most  desirable  that  the  attacking  aircraft  approach 
the  target  along  the  same  line  as  the  illuminator  in  three-dimensional 
space.  The  reflective  cross  section  will  be  higher  because  the  diffuse 
reflective  lobe  of  the  surface  material  will  be  greater,  and  the  receiver 
will  be  viewing  more  illuminated  surface. 

Another  effect  is  the  reflectivity  of  the  background.  Grass,  for 
example,  exposes  a  much  larger  reflect’’  «.g  surface  when  viewed  and 
illuminated  from  a  shallow  angle.  The  1.065-[i  bi-directional  reflectance 
of  grass  varies  from  50%  (receiver  overhead)  to  80%  (receiver  at  a  10° 
dive  angle)  when  illuminated  at  a  27°  dive  angle  (Kef  12:909).  Roads, 
especially  wot,  have  a  much  more  specular  reflection  at  low  angles  of 
illumination  and  show  less  backscatter  than  the  surrounding  fields.  If 
the  local  background  lias  a  much  higher  reflectance  than  the  target  and 
is  within  an  acceptable  kill  radius  of  the  target,  it  will  be  better  to 
illuminate  the  background  material  beside  the  target. 

Seeker  Receiver  Characteristics  (In  General) 

Laser  seeker  system.'--  are  classified  by  that  characteristic  which 
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most  limits  their  usefullness  to  the  target-acquisition  problem.  If  a 
seeker  is  not  cryogenically  cooled,  its  detector  housing  or  load 
resistor  thermal  noise  will  be  the  limiting  factor  (internal-noise- 
limited).  If  the  detector  has  a  low  efficiency  <?.c  the  wavelength  used 
or  is  followed  by  low-gain  amplifiers,  then  it  is  sensitivity-limited. 
Either  of  these  are  termed  seeker-limited  systems. 

The  newest  cryogenically  cooled  (4°K)  avalanche  photodiodes  and 
photomultiplier  tubes  approach  the  theoretically  possible  conversion 
efficiencies;  but  then  these  are  background-limited  by  the  flux  from 
the  terrain  within  the  seeker  field  of  view  (5-500  milliradians).  The 
radiance  of  the  background  may  be  several  orders  of  magnitude  less  than 
the  target,  but  due  to  the  huge  viewed  area,  may  give  more  net  flux  at 
the  detector.  This  can  be  improved  by  either  using  expensive,  high- 
resolution  imaging  systems  or  narrowing  the  field  of  view  until  the 
statistical  probability  of  the  target  not  being  in  the  field  of  view 
increases  unacceptably.  This  constant  solar  illumination  results 
primarily  in  a  largo  DC  component  of  the  detector  signal,  which  can  be 
filtered  out.  The  white  noise  component  of  the  signal,  however, 
requires  considerable  filtering  and  electronic  analysis  to  detect  the 
desired  signal. 

Most  detector  heads  currently  used  arc  non-cooled,  so  are  seeker- 
limited  systems.  The  signal/noise  ratio  cf  these  detectors  varies 
considerably  with  the  temperature  of  the  detector  and  its  load  resistor. 
An  aircraft  from  a  hot  tropical  base  for  a  low  altitude  flight  to  a 
local  target  will  have  a  detector  system  at  100°-150°F,  depending  on  how 
long  the  weapons  have  been  sitting  in  the  sun.  The  same  system  cruising 
for  two  hours  at  AO, 000  feet  will  have  a  detector  cold-soaked  to  -50°F. 
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Unless  compensated  internally  for  this  temperature  change,  the  200°F 
cooler  temperature  will  have  a  beneficial  effect  of  several  orders  of 
magnitude  on  the  signal/noise  ratio  of  the  seeke:  and  the  sensitivity  of 
the  system,  and  hence  on  the  lock-on  range.  Internal  temperature 
compensation  is  therefore  necessary  for  prediction  studies  and  analysis. 

Whether  background-limited  or  seeker-limited  the  system  can  be  as 

1.  contrast-limited  problem  as  in  the  electro-optical  target  seekers 
which  use  a  high-resolution  TV  imaging  system,  or 

2.  "one  brightest  spot  in  the  field  of  view"  problem,  or 

3.  "centroid  of  entire  illumination  in  the  field  of  view"  problem. 

System  1  is  the  most  electronically  complex  and  expensive,  and  has  the 
disadvantage  thac  the  operator  must  designate  the  target  with  a  pointing 
device  in  the  same  imaging  system.  Systems  2  and  3  have  the  problem 
that  the  portion  of  the  laser  energy  falling  into  the  more  reflective 
background  may  create  a  brighter  return  than  the  desired  target.  The 
same  problem  is  created  when  haze,  clouds,  or  foliage  is  in  the  beam 
between  the  illuminator  and  the  target.  The  seeker  will  either  track  the 
brightest  of  the  targets  available  or  the  centroid.  As  a  centroid 
seeker  gets  closer  to  the  target,  it  will  eventually  reach  a  point  whore 
one  of  the  bright  spots  reaches  the  edge  of  its  field  of  view,  and  then 
will  change  course  to  attempt  to  track  the  remaining  target  or  targets. 

This  problem  can  bo  considerably  alleviated  by  designing  arrays  and 
logic  circuits  that  discriminate  against  extended  returns  from  clouds, 
haze,  and  flat  terrain,  and  discriminate  for  first/last  pulse  of  the 
proper  strength  and  shape  to  bo  from  a  sudden,  hard  surface  too  large  to 
be  tree  branches  and  too  small  to  be  a  sloping  hillside. 

Seeker-limited  systems  are  limited  by  their  lack  of  sensitivity, 
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and  high  internal  noise;  and  require  illuminances  on  the  order  of  10 

-5  2 

to  10  watts/cm  .  Background-limited  systems  can  work  with  illuminances 
on  the  order  of  10  ^  to  10  10  watts /cm4-  with  comr.ercially  available 
cryogenically  cooled  detectors.  The  latter  are  more  expensive,  require 
considerable  field  maintenance,  and  are  not  rugged  for  field  use,  but 
if  the  problems  with  pulse-shape-analysis  and  multi-target  discrimination 
can  be  solved,  they  can  extend  the  useful  seeker  clear-air  lock-on  range 
by  the  square  root  of  the  increase  in  sensitivity.  In  this  author’s 
opinion,  this  is  the  present  serious  limiting  factor  in  extending  the 
lock-on  ranges  of  laser-guided  systems. 

This  chapter  has  discussed  systems  in  general  and  not  any  one 
operational  system  specifically.  The  information  is  a  general  summary 
of  the  characteristics  and  problems  of  system  types  in  general,  based 
on  over  300  unclassified  reports  and  articles.  This  completes  the 
chapter  on  the  geometry  and  design  of  laser  systems  in  the  absence  of 
an  attenuating  atmosphere.  The  next  chapter  will  deal  with  the  physical 
nature  of  that  aspect  of  the  lock-on  range  problem. 
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III.  Attenuation  Mechanisms  &  Atmospheric  Structure 


Attenuation  Mechanisms 

Laser  beams  in  a  real  atmosphere  are  attenuated  by  absorption  (both 
band  &  line,  in  both  gases  &  aerosols  in  the  optical  path)  and  scattering 
(Rayleigh,  Mie,  and  turbulence  cells).  This  chapter  will  discuss  these 
mechanisms  and  their  effects  on  1.065-p,  laser  beams  in  the  bottom  10  Km 
of  the  atmosphere. 

Attenuation  by  Gases.  Water  is  the  main  absorbing  gas  at  1.065  p 
in  the  lower  atmosphere.  Only  in  special  industrial  areas  with  a  strong 
thermal  inversion  may  there  be  enough  concentration  of  specialized  waste 
gases  to  cause  significant  attenuation  at  certain  wavelengths.  Rayleigh 
scattering  at  1.065  (i.  is  negligible  (Ref  37:7-30). 

From  McClatchey  (Ref  26: 36, 59, A 3)  and  Eldridge  (Ref  17:710)  one  cm 
of  precipitable  water  vapor  will  have  a  transmission  of  96%  +  1/2%  at 
1.065  p.  Since  the  attenuation  coefficient  B  is  defined  by: 

pa¥106io(7)  (2*> 


where  R  is  the  path  length  and  T  is  the  transmission,  this  yields  an 
attenuation  coefficient  Q  =  0.175  +  0.022  dB/cm  of  prc-cipitable 
water  vapor  in  the  optical  path.  The  water  content  of  a  path  can  be 
obtained  from  the  relative  humidity  and  temperature  by  Figure  8  or  as 
follows.  By  definition, 


•  .  G 

Relative  Humidity  =  f  =  — 

e 

s 


existing  water  vapor  pressure 
saturated  water  vapor  pressure 
for  that  temperature. 


(25) 


From  Haurwitz  (Ref  21:9): 


e  =  (A. 58)  10 

o 


^t  +  2iy 


(26) 


where  t  is  the  temperature  in  degrees  C  and  es  =  saturated  vapor  pressure 
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Fi(j.  8.  Watei  Vapor  Concentration  per  Km  of  path  as  a  Function  of  Temperature 
and  helative  Humidity  (Kef  26:A5) 
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in  nun  Hg.  This  would  match  the  meteorological  data  (Ref  26:A5)(Ref  40: 

49)  better  if  the  exponent  were  8.3  instead  of  7.5j  therefore  the 

numeric  8.3  will  be  used  in  this  thesis.  From  Byers  (Ref  7:152,154), 

the  absolute  humidity  (AH)  is  given  by: 
em 

iii  _  w  _  <6--  e 

R(t  +  273)  "  R  (t  +  273) 

where  m  is  the  molecular  weight  and  R  is  the  universal  gas  constant, 
w 

The  number  of  cm  of  precipitable  water  in  a  column  of  air  is  the 
absolute  humidity/10  (Ref  26:A3),  so  combining  the  above  information: 


_  .622  e _  _  .285  f  10 

CmH20  10  K(t  +  273)  R(t  +  273) 


/8«3  t  v  .80  3  t  v 

.285  f  10  t  +  2/3  _  (136)f 10  +  2/3’ 


t  +  273 


The  attenuation  coefficient  at  1.065  ji.  due  to  water  vapor  is  therefore: 

=  (.175)  (cm^)  in  dB/Km  (29 


Attenuation  by  Turbulence  Cells.  Turbulence  cells  (small  pockets 
of  hotter  or  cooler  air  rising  or  settling  turbulently  due  to  thermal 
convection  in  the  atmosphere)  by  themselves,  will  not  absorb  laser  beam 
radiation,  but  they  will  scatter  and  defocus  the  beam  locally.  From  one 
point  in  the  beam  at  the  target,  the  laser  source  will  appear  to 
scintillate  about  the  average  non-turbulent  intensity.  There  will  be  some 
beamspreading,  which  was  found  by  Buck  (Ref  5:705)  to  be  on  the  order  of 
d  =  aR  ’  for  \  experiments  at  Boulder,  Colorado,  using  an  a  =  4.5 
microradian  beam  divergence  over  ranges  from  0.5  Km  to  145  Km  in 
relatively  clear  air,  in  mountains  2  Km  above  sea  level.  Under  hazy 
conditions  and  daylight  hours  the  exponent  of  R  increased  (Ref  5:705). 

Turbulence  cells  will  cause  a  statistical  variation  in  the  beam 
intensity,  yielding  a  pattern  analogous  to  Figure  2.  The  large  target 
area  will  average  those  out,  so  ’-he  total  energy  reflected  by  the  target 
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will  not  fluctuate  significantly.  The  small  seeker  detector  area, 
however,  will  not  average  the  fluctuations  caused  by  the  turbulence 
between  the  seeker  and  the  target,  so  there  will  be  statistical 
fluctuations,  up  to  200  Hz  (Ref  5:706)  in  the  receiver  signal. 

Due  to  the  stated  assumptions  in  pointing  &  tracking,  beam  wander  is 
disregarded. 

Attenuation  by  Aerosols.  Absorption  and  scattering  depend  strongly 
on  the  aerosol  composition  and  particle  size  distribution.  The  oceanic 
salts  are  the  primary  constituents  of  the  soluble  condensation  nuclei, 
containing  primarily  the  Ca,  K,  C03,  Br,  Cl,  Sr,  and  Na  compounds.  To 
this  is  added  huge  amounts  of  insoluble  materials  from  land  soils, 
aromatic  and  organic  materials  from  vegetation,  and  products  of  a  very 
complex  photochemical  reaction  of  all  these.  The  mixing  ratios  of  these 
compounds  vary  considerably,  as  there  is  rapid  transport  of  aerosols 
for  at  least  1000  miles  (Ref  24:  sect  2.3)  Apparently  there  are  no 
purely  continental  or  purely  maritime  aerosols,  but  growth  curves  of 
aerosol  particles  with  increasing  relative  humidity  indicates  that  the 
condensation  nuclei  for  central  Europe  average  roughly  70  -  80^ 
insoluble  materials  (Ref  24:169).  The  condensation  nuclei  have 
specialized  composition  around  industrial  regions,  heavily  vegetated 
regions,  mountains,  dusty  areas,  and  windy  areas.  It  takes  several  days 
for  sufficient  mixing  to  smooth  out  the  composition  mixing  ratios. 

Attenuation  by  solid,  dry  aerosols  is  a  mixture  of  absorption  and 
scattering.  CRL  aerosol  analysis  data  (Ref  26:9,21)  computes  the 
attenuation  coefficient  for  1.065  a  at  23"  absorption  and  77%  scattering 
for  visibilities  above  5  Km  and  relative  humidities  below  70". 
Attenuation  by  wot  aerosols  (aqueous  solutions  or  a  thick  layer  of  water 
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around  an  insolublo  solid  coro)  is  almost  ontiroly  scattering.  Above 
70%  relative  humidity  the  wet  aerosol  particles  rapidly  dominate.  During 
decreasing  humidities  the  supersaturated  aqueous  solutions  remain  liquid 
down  to  about  40%  relative  humidity. 

The  index  of  refraction  of  water  at  1.065  p  is  1.325  (Ref  10:244), 
and  of  a  large  proportion  of  materials  found  in  aerosol  particles  is 
between  1.3  and  1.5  (Ref  13:1)  with  some  between  1.5  and  1.7  (Ref  41: 
328).  Therefore  the  real  index  of  refraction  of  aerosol  particles  will 
be  between  the  extremes,  and  will  bo  well  towards  the  lower  value  for 
the  higher  humidities,  where  the  proportion  of  water  in  the  aerosol  will 
be  very  large  (Ref  29:6).  This  has  been  validated  by  attenuation 
measurements  and  calculations  by  Gibbons  (Ref  19:172-176)  within  +  4% 
of  the  theoretical  values  for  most  nonindustrial  hazes. 


Wet  aerosols  are  spherical;  the  highly  irregular  shape  of  the  dry 
aerosols  will  average  out  to  a  mean  effective  spherical  radius  since  any 
possible  orienting  influence  of  the  pulsed  laser  beam  is  randomized  by 
the  Brownian  thermal  motion  around  the  aerosols.  This  is  consistent  with 
the  dry  dust  scattering  properties  measured  by  Volz  and  Goody  (Ref  38: 
390).  The  absorption  coefficient  for  spherical  aerosol  particles  is 


approximately  (Ref  29:20,22): 

r _ 


P  ,  ='TT' 

aerosol 


absorotion 


il  r  1  +  +  —  r — n(r  )  dr 

3  P  2xi1’  ftv2  ,2  P 


whore : 

n(i’p)  =  particle  size  distribution  function, 
r  =  radius  of  aerosol  particle 

P  -r 

n*  =  imaginary  index  of  refraction  of  the  aerosol,  between  10~°  and  10 
for  most  water- solution  aerosols. 

2  Tf  T 

x  =  — — p  and  X  =  wavelength  of  laser  radiation  =  1.065  p. 


The  4/3  in  Eq  (30)  adjusts  the  a,  pLcximation  for  a  better  fit;  of  the 
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rigorous  and  exact  Mie  solutions.  The  scattering  coefficient  for 
spherical  aerosol  particles  is  given  by  (Ref  29:25;: 


8  _  n  =  IT  f  r2  n(r  )  Q  (m,r  ,\)  dr 

Htotal  J  p  p  xscatt  p  p 

scattering  min 


where : 

n(r  )  =  particle  size  distribution  function. 

0  ^(m.r  ,X)  =  Mie  scattering  efficiency  factor. 

^scatt  p 

m  =  real  index  of  refraction  of  the  aerosol. 

Many-particle  scattering  can  be  neglected. (Ref  29:23)  because  even 
in  the  densest  recorded  fogs,  the  particles  are  still  separated  by 
distances  over  20  times  their  radius.  From  his  measurements  Zuev  (Ref  42: 
140)  determined  that  below  9  optical  depths  (1  optical  depth:  length  at 
which  intensity  decreases  to  1/e  due  to  absorption  &  scattering),  single 
scattering  is  valid  for  narrow  collimated  laser  beams.  Above  9  optical 
depths  (=39  dR  attenuation  =  1/SOOO  origionul  intensity)  multiple 
scattering  may  be  significant  for  targets  more  than  9  times  the  diameter 
of  the  laser  beam. 

Additionally,  the  energy  in  the  forward-scatter  lobe  broadens  the 

beam  and  decreases  8  .  ♦  For  small  particles  (haze  and 

toucvi  scacccxmg 

most  clouds)  the  forward-scatter  lobe  is  extremely  broad,  and  hence  the 
energy  scattered  onto  the  distant  target  is  small.  Snow  also  has  a  very 
broad  lobe.  Large  rain  particles  have  a  strong  and  narrow  forward- 
scatter  lobe  that  decreases  the  total  attenuation  coefficient  in  rain  by 
8.4%  at  0.6323  p  and  32.6%  at  3.5  pi  (Ref  11:733).  Since  he  scattering 
raindrops  were  much  larger  than  the  wavelengths  cited,  the  scattering 
efficiency  factor  Q  is  constant  with  wavelength,  and  therefore  a  linear 
interpolation  is  valid.  Thu?-  aerosol  scattering  is  decreased  by  13%  at 
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Table  11 


Beamspreading  as  a  Function  of  Attenuation  and  Type  of  Weather 


Weather  Precip. 

Rate 


Beam  Total- 
Divergence  Attenu; 
(milliradians)  dB/Km 


.209  +  .005 
.295  +  .005 

.48  +  .02 

(extrapolated) 


Total-beam  Effective  Attenuation 
Attenuation  at  the  Beam  Center 
dB/Km  Caused  by  Broadening 


1.1  dB/Km 
2.8  dB/Km 


0.85  dB/Km 
1.6  dB/Km 


(From  Ref  11:753-4) 


Fig  9.  Bcamspi ending  as  a  junction  o»  Kent  her  Attenuation 
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1.065  p.  in  rain. 

Chu  &  Hogg  (Eef  11:753-4)  measured  beam  broadening  of  a  gaussian 
beam  over  a  2.6  Km  path  at  0.6328  p  in  rain,  fog,  and  snow  for  different 
attenuation  rates.  The  unattenuated  beam  had  a  whole-angle  divergence 
of  0.2  milliradians.  The  results  are  in  Table  II  and  Figure  9.  In 
general,  the  larger  the  beam  divergence,  the  larger  the  effective 
forward  scattering.  This  effect  could  be  used  effectively  when 
illuminating  large  surface  areas  through  rainstorms. 

Particle  Size  Distributions 

To  solve  Eqs  (30)  &  (31)  it  is  necessary  to  know  the  particle  size 

distribution  n(r  ),  and  its  maximum  and  minimum  r  and  r  .  . 

p  max  min 

Unfortunately  our  knowledge  in  this  field  is  very  primitive.  The 
particle  size  distribution  has  been  traditionally  determined  by  means 
that  inherently  disrupt  the  distribution  itself  or  the  gathering 
efficiency  changes  l'elativo  to  particle  size  in  an  unknown  manner. 

The  first  method  was  to  mount  a  flat,  cold,  dry  plate  and  a  flat, 
cold,  oily  or  sticky  plate  (or  box)  outside  an  aircraft,  facing  the 
airflow,  and  fly  through  the  haze  or  cloud.  The  collected  aerosols  wore 
viewed  under  an  optical  or  an  electron  microscope.  Unfortunately  the 
very  small  particles,  under  one  micron  radius  and  therefore  exactly  in 
the  realm  of  interest,  would  follow  the  airflow  around  the  plate  or  box, 
resulting  in  a  smooth  size  distribution  favoring  the  heavy  particles. 
Secondly,  the  particles  would  evaporate  in  the  electron  microscope 
chamber  vacuum  or  under  the  high  intensity  light  or  electron  beam. 

The  second  method  was  to  view  the  aerosol  particles  directly 
through  a  microscope  mounted  at  the  neck  of  the  venturi  tube  on  the  side 
of  an  aircraft.  Unfortunately  t*”'  considerable  change  in  temperature, 
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pressure,  and  spacr.ng  of  the  aerosol  particles  caused  rapid  evaporation 
of  the  smaller  liquid  particles  and  coagulation  of  the  larger  particles. 

The  third  met!  od,  determination  of  size  distributions  by  measure¬ 
ments  of  attenuation  of  several  laser  beams  through  the  same  path, 
requires  a  prior  assumption  of  the  general  form  of  the  distribution 
function.  A  variation  of  this  involves  measuring  the  scattered  light 
from  the  beam  at  different  angles  to  determine  a  total  scattering 
function  for  all  the  particles.  Both  yield  generally  acceptable 
approximations,  but  require  prior  assumptions  about  the  unknown  size 
distribution  being  tested.  A  method  needs  to  be  developed  to  directly 
study  the  aerosols  at  a  distance  from  the  aircraft  that  will  yield  a 
direct  measurement  of  undisturbed  aerosols. 

Particles  under  0.1  p  are  inefficient  scatterers.  Particles  over  a 

20  p  radius  gradually  precipitate,  with  settling  rates  proportional  to 

3 

mass  (proportional  to  r  )  (kef  24:111).  Fog  droplets  range  between  1  and 
100  p,  while  raindrops  average  1000  p.  For  the  haze  particles  between 
0.1  and  100  p,  the  expression  (Ref  29:10): 


n(y 2 


dN 

dr 


k 


where: 

N  =  number  of  particles  with  radius  r  , 

P 

k,b,6,y  are  experimentally  determined  constants, 


(32) 


is  generally  agreed  upon  as  representat.ing  a  final,  steady-state  size 
distribution  and  a  first  approximation  to  most  real,  atmospheric  aerosols. 
Unfortunately  the  local  condensation  nuclei  input  and  the  dynamic 
atmosphere  prevent  a  steady  state  condition  from  existing.  Bimodal 
(double-peaked)  distributions  are  relatively  common  (Ref  24:140)  (Ref  17: 
18).  It  may  take  several  days  for  particles  under  1  p  to  reach  single- 
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mode  equilibrium  (Ref  24:130).  Chu  & 

a  =  r  /r  ,  where  r  is  the  radius  of 
pm  n 

number  density.  Eq  (32)  then  is  (Ref 


Hocg  modified  Eq  (32)  by  defining 
the  particle  size  with  the  highest 
11:730) 


n(a)  =  k  a6  e~l'aY 


(33) 


This  representation  is  more  flexible  and  in  agreement  with  measurements 
from  a  wide  range  of  meteorological  conditions  (Ref  29:12). 

Below  99 %  relative  humidity,  the  r^  is  less  than  0.1  p,  and  the 
decrease  in  number  density  with  increase  in  radius  can  be  approximated 
by  (Ref  24:116,124): 

1M  —  ~  A\ T  —  A 

(34) 


dN _ rp  dN  _  -g 

d(log  r')=  drp  crp 


dN  -d 

or  —  =  c  r 
dr  p 

P 


where  c,  d,  and  g  are  empirical  constants. 

Above  99%  relative  humidity,  the  condensation  process  overwhelms 
the  evaporation  process  for  small  particles,  and  the  coagulation  process 
dominates  for  the  larger  particles,  so  the  rm  increases  rapidly  from  . 
below  0.1  p  to  any  value  between  0.1  and  1000  p.  Rain  and  fog  readily 
form.  For  heavy  drizzle  and  rain,  a  gaussian  distribution  is  frequently 
used  (Ref  29:15).  When  a  large  majority  of  spherical  particles  is  over 
5-10  times  the  laser  wavelength,  a  total  scattering  efficiency  Q  of 
about  2  can  bo  assumed  in  the  Mie  scattering  equations  for  wide 
distribution  patterns,  thus  greatly  simplifying  numerical  calculations. 
Measured  and  calculated  particle  size  distributions  are  given  in 
References  17,  23,  29,  and  others. 

It  must  bo  emphasized  that  the  real-world  atmospheres  are  very 
complex,  and  the  study  of  weather  and  particle  size  distributions  is 
still  in  its  infancy.  Attenuation  predictions  through  aerosols  from 
meteorological  measurables  can  be  only  crudely  approximated. 
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i. 

Atmospheric  Models  (Vertical  Structure) 

The  particle  size  distribution  does  not  remain  constant  as  altitude- 
increases*  The  upper  limit  of  suspended  aerosol  size  decreases  from 
about  20  (i  at  sea  level  to  10  u  at  3  Km  altitude  (Ref  24:118),  with  a 
corresponding  decrease  in  the  particle  distribution  between  2  u  and  the 
upper  limit  for  that  altitude.  Below  2  pi  the  distribution  remains  about 
the  same.  This  settling-out  is  particularly  noticeable  for  an  atmosphere 
with  no  thermal  mixing  and  no  thermal  inversions  for  several  days.  The 
heavier  aerosols  will  settle  down  to  an  altitude  where  the  increased 
density  (more  Brownian  motion)  will  keep  them  aloft.  The  aerosol 
profile  will  eventually  conform  to  the  CRL  exponential  standard 
atmosphere. 

If  the  bottom  region  of  the  troposphere  up  to  the  first  strong 
thermal  inversion  has  considerable  convective  or  thermal  mixing,  the 
sedimentation  process  of  large  particles  will  be  negated,  and  the  result 
is  a  homogeneous  atmosphere.  Indeed,  several  cases  were  measured  by  US 
aircraft  making  whole-volume  scattering  measurements  over  Germany  (Ref 
18),  showing  a  constant  scattering  coefficient  up  to  the  strong 
inversion  layer,  between  1  and  3  Km  altitude  on  different  days,  and  then 
a  sudden  decrease  over  the  range  of  200  meters  in  the  inversion  to 
values  an  order  of  magnitude'  lower  than  below  the  inversion.  The 
measurements  above  the  inversion  showed  either  an  exponential  decreas . 
characteristic  of  a  steady-state,  no-mixing  atmosphere,  or  a  nearly 
constant  and  low  value  above  that  inversion  layer. 

Frequently  in  the  early  or  mid-afternoon,  the  first  strong 
inversion  layer  disappears  due  to  thermal  heating  from  below,  and  the 
formerly  trapped  aerosols  dissipe'-o  rapidly  upwards  to  the  next 
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inversion  layer.  Equilibrium  is  then  established  within  a  few  hours. 
This  process  may  bo  progressively  repeated  through  successively  higher 
inversions.  The  atmosphere  during  the  transition  is  somewhere  between 
the  above  mentioned  two  extreme  atmospheric  states. 

For  calculations  in  this  thesis,  it  will  be  assumed  that: 

1.  The  atmosphere  is  the  same  regardless  of  azimuth  from  the  target, 
for  10  Km  radius. 

2.  The  atmosphere  may  be  horizontally  stratified,  exponential  and 
unstratified,  or  homogeneous.  The  first  case  will  be  modeled  after 
flight  test  data  given  by  Dr.  Robert  W.  Fenn  (Ref  IS),  from  his  analysis 
of  data  of  Flight  C139,  Haven  View,  on  3  June  1970,  near  Memmingen, 
Germany.  The  second  case  will  be  modeled  after  AFCRL  field  measurements 
(Ref  37:7-1).  The  third  case  will  be  modeled  after  the  results  of  the 
first  case. 

3.  The  atmosphere  everywhere  within  the  10  Km  radius  of  the  target  has 
a  relative  humidity  below  99%:  no  fog  or  stratus  layers  exist  with  a 
different  particle  size  distribution. 

The  first  condition  is  very  restrictive.  The  atmosphere  in  any  one 
layer  is  inhomogeneous.  The  calculations  will  give  an  average,  best, 
or  worst  case  depending  on  the  corresponding  attenuation  coefficient 
used.  Multi-transmissomcter  studios  at  airfields  at  Atlantic  City,  New 
York  City,  Chicago,  Denver,  and  Los  Angeles  (Ref  31)  show  that  when 
visibilities  are  reduced  to  2  Km  or  less,  that  25%  variations  often 
occur  in  atmospheric  ti’ansmissivitics  between  successive  readings  15 
seconds  apart  at  the  same  point,  even  when  the  transmissometer  averages 
over  15  seconds  observation  time.  Furthermore,  the  three  transmissom- 
eters  along  a  3  Km  runway  path  ,'how  incoherent  fluctuations,  as  the 
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inhomogeneties  for.n,  dissipate,  and  drift  about.  When  the  visibility  is 
below  0.5  Km,  6002  fluctuations  often  occur.  Above  3  Km  visibility, 
the  local  variations  average  out  because  the  observer  can  see  over  a 
longer  path  and  the  attenuation  coefficients  are  much  lower.  These 
tests  were  made  with  light  in  the  visual  range,  but  the  results  are 
equally  applicable  to  the  very  nearby  1.065  p  interest  of  this  thesis. 

The  last  condition  is  also  very  restrictive.  As  the  aerosols  rise 
above  the  inversion  layer,  the  temperature  drops  rapidly  and  the 
relative  humidity  increases  suddenly.  Thus  stratus  layers  frequently 
form  at  the  top  of  the  inversion  layer,  wich  resulting  severe 
attenuation  for  those  few  hundred  meters  in  the  stratus  layers. 

Even  with  these  three  restrictions,  useful  and  meaningful 
information  can  be  obtained  from  the  evaluation  of  the  three  cases  of 
atmospheric  models  listed  above.  With  these  models  of  the  vertical 
profile  of  the  atmosphere  and  the  models  of  the  particle  size 
distribution,  it  is  now  possible  to  combine  them  and  determine  the 
attenuation  of  a  laser  beam  in  a  real  atmosphere,  if  all  the  previously 
mentioned  empirical  numbers  are  known  as  a  function  of  measurable 
meteorological  parameters. 

Attenuation  Coefficients  and  Weather  Phenomena. 

At  1.065  ps 


(From  Eq  (30) 


^total  attenuation  =  ^  ^li^O 

+  B 

^aerosol 
absorption 

**"  ^aerosol  (rrom  Eq  (31) 
scattering 


(35) 


Since  raindrops  average  1000  p  in  radius,  much  larger  than  1.065  p,  the 
scattering  attenuation  (corrected  for  forward  scattering)  is  (Ref  22:6-1 


8  ,  =  0.88  h 

•aerosol  scattering 


0.683 


dB/Km 


(36] 
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.where  h  =  mm  of  rain  per  hour*  Since  rain  is  almost  totally  pure  water, 

B  ,  ,  .  is  much  less  than  {3  ,,  ,  , .  .  .  . 

aerosol  absorption  aerosol  scattering 

The  attenuation  properties  of  snow  (almost  completely  scattering) 
are  complex,  but  in  general  the  attenuation  for  dry  powder  snow  is  much 
less  than  that  of  an  equivalent  amount  of  water  in  the  form  of  rain.  Wet 
snow,  however,  has  a  thin  liquid  water  surface  over  an  irregular  surface 
much  larger  than  the  laser  wavelength,  attenuates  much  more  than  an 
equivalent  amount  of  rain,  and  falls  more  slowly  (Ref  22:6-14)* 

The  total  atmospheric  attenuation  is  the  integral  of  Eq  (35)  along 
the  entire  optical  path.  A  correlation  of  n(r^)  with  something  readily 
measurable  by  weather  observers  is  not  yet  known;  but  within  certain 
limits  of  weather  and  wavelength  some  approximations  can  be  made. 

For  haze  with  relative  humidity  less  than  9 9%,  the  particle  size 
distribution  can  be  crudely  classified  as  continental  or  maritime  or  a 
combination  of  the  two,  as  defined  by  Junge  and  modified  by  Chu  &  Hogg. 
Figure  10  plots  the  computed  total  aerosol  attenuation  coefficients  from 
the  CRC  Handbook  of  Lasers  (Ref  40:43),  CRL  Optical  Properties  of  the 
Atmosphere  (Ref  26:15,16,19,21)  and  CRL  Vertical  Attenuation  Model  (Ref 
33:28-47).  All  are  representative  of  typical  aerosols. 

If  aerosol  scattering  is  the  primary  attenuating  mechanism,  perhaps 
a  way  around  the  lack  of  knowledge  of  precise  particle  size  distribu¬ 
tions  is  to  measure  the  attenuation  at  two  wavelengths  relatively  close 
(compared  to  the  breadth  of  the  distribution)  to  the  1.065  p  laser 
wavelength.  Since  both  the  continental  and  maritime  hazes  (and  therefore 
a  mixture  of  the  two)  cover  the  range  from  0.02  -  10  p  (and  fogs  up  to 
100  p)  a  reasonable  choice  for  the  other  wavelengths  would  be  at  0.42  p 
and  0.695  p  in  the  visible  spec*  *un.  Measurements  of  attenuation  at 
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Fie*  10*  Computed  Aerosol  Attenuation  Co-.-M icieuc-  iron  Three  Sourc 
(Ref  40:43)  (Rof  33i 2b-^i7>  (Ref  2o:15-2l) 
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these  wavelengths  can  be  measured  by  measuring  t.ae  meteorological  range 
through  narrowband  filters  at  these  wavelengths.  If  the  particle  size 
distribution  is  broad  compared  to  the  spacing  between  0.42  and  1.065  p., 
does  not  have  sharp  peaks  near  or  between  these  wavelengths,  and  is  not 
multi-modal  (or  if  so,  it  has  all  modes  above  one  micron),  then  there 
should  be  some  simple  relationship  between  the  attenuation  at  the  three 
wavelengths  above.  This  criteria  are  met  by  almost  all  hazes,  fogs, 
and  clouds. 

Meteorological  range  (visibility,  or  meteorological  visibility)  in 
white  light  is  the  distance  at  which  the  observer  can  barely  detect  a 
very  large,  black  object  against  a  white  background  (apparent  contrast 
is  reduced  from  unity  to  2%).  Due  to  a  lack  of  information  about  the 
contrast  threshold  of  the  photopic  eye  as  a  function  of  wavelength,  it 
will  be  assumed  that  the  eye  can  detect  this  21  contrast  in  narrowband 
red  light  and  narrowband  blue  light  as  well  as  in  white  light.  This 
yields  a  relation  between  the  attenuation  coefficient  £  and  the 
meteorological  visibility  V  (from  Appendix  D): 

3  =  in  dB/Km 

true 

The  reported  visibility  from  weather  observers  is  generally  about 
75^  of  the  true  meteorological  range  (Ref  18)  because! 

1.  Perfectly  black  objects  on  a  perfectly  white  background  are  ideal 
cases,  nonexistant  in  field  conditions. 

2.  There  are  a  limited  number  of  objects  around  the  observer,  statisti¬ 
cally  reducing  the  chances  of  one  being  exactly  at  the  meteorological 
range.  The  next  object  inside  the  meteorological  range  will  be 
reported  if  its  inherent  contrast  is  close  enough  to  unity. 

3.  The  observer  reports  the  objects  he  can  identify.  The  contrast 
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threshold  for  recognition  is  higher  (about  5%)  than  for  detection 
(defined  at  2%).  Therefore  Eq  (37)  becomes: 


visually  measured 

where  Vyisually  measured  =  reported  visibility  =  apparent  visibility. 

''  The  reported  (rather  than  true)  visibility  will  be  used  in  this  thesis. 

If  the  true  visibility  were  to  be  used,  the  spectral  characteristics  and 
inherent  contrast  of  all  visible  objects  around  the  observer  for  all 
weather  conditions  and  sun  angles  must  be  tabulated,  so  when  the 
observer  reports  that  he  can  see  object  X,  the  true  meteorological  range 
could  be  determined  from  the  previously  tabulated  data.  This  is  never 
done,  and  it  is  the  intent  of  this  thesis  to  arrive  at  a  methodology  by 
which  an  observer  without  considerable  equipment  can  provide  information 
which  will  allow  the  determination  of  the  attenuation  coefficient  of 
1.065  (i  laser  "beams  in  real  atmospheres. 

The  irradiance  from  the  sun  is  not  uniform  with  wavelength,  but 
decreases  by  30%  at  0.42  and  by  20%  at  0.695  {t,  compared  to  peak  yellow 
irradiance  (Ref  37:16-2).  The  sensitivity  of  che  eye  at  0.42  a  and  0.695 
[1  is  less  than  1%  of  the  peak  yellow-green  value  (Ref  20:84l)-i  This 
severe  decrease  in  perceived  light  will  make  meteorological  range 
determinations  more  difficult  and  may  change  the  contrast  threshold  of 
the  eye. 


Rayleigh  scattering  must  be  subtracted  from  the  attenuation  observed 
at  blue  wavelengths  and  water  vapor  absorption  added  to  the  attenuation 


at  1.065  [i  to  make  the  attenuation  coefficient  comparison  meaningful. 

At  0.42  p,  Rayleigh  scattering  is  0.16  dB/Km  at  sea  level  (Ref  16:7-5). 
Water  vapor  absorption  will  range  from  near  zero  to  0.9  dB/Km. 
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Attenuation  coefficients  computed  by  Carrier  (Ref  8:1212)  as  a 
function  of  wavelength  for  clouds  show  a  series  of  smooth,  flat  curves 
from  0.488  to  1.065  a,  plotted  in  Figure  11.  The  size  distribution 
functions  of  lines  i,  2,  and  3  have  small  but  noticeable  tails  of  large 
droplets  over  6  p;  the  functions  for  lines  4,  5,  6,  and  7  have  smaller 
and  less  significant  tails;  and  line  8  has  practically  no  tail.  Note  the 
difference  that  the  tail  of  large  droplets  makes  in  the  total  attenuation 
coefficient  for  all  wavelengths.  Almost  all  the  particles  are  between 
1  and  15  p  in  radius.  Rayleigh  scattering  and  water  vapor  absorption 
are  negligible. 

Attenuation  coefficients  were  computed  by  Rensch  (Ref  29:39,40)  for 
a  continental  model  aerosol  with  relative  humidities  below  90%  or  at  95%. 
The  computations  were  made  for  several  slopes  of  the  distribution  curve 
from  0.1  -  1  p;  the  results  are  plotted  in  Figure  12.  For  relative 
humidities  below  90%  the  attenuation  at  1.065  p  can  be  determined  by  a 
linear  graphical  extension  of  the  attenuation  coefficients  at  0.42  p  and 
0.695  p  when  the  0.42  p  readings  are  corrected  for  Rayleigh  scattering 
and  water  vapor  absorption  is  added  for  the  1.065  p  attenuation 
coefficient.  At  95%  relative  humidity,  a  linear  extrapolation  of  the 
attenuation  coefficients  results  in  an  overestimation  of  the  attenuation 
coefficient  at  1.065  p  by  approximately  0.35  dB/Km. 

Similar  computations  by  Blattner  (Rtf  4:29)  using  the  Elterman  haze 
model  for  meteorological  ranges  from  3  Km  -  50  Km  show  that  a  linear 
extrapolation  from  0.45  p  through  0.65  p  to  1.065  p  is  reliable.  Other 
calculations  by  Doirmondjian  (Ref  14)  and  Cato  (Ref  9)  concur  for 
relative  humidities  below  90%  or  above  saturation  (clouds  and  fogs).  It 
is  in  the  transition  range  between  90  -  100%  that  th.  aerosol  particle:; 
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. 3  .5  .7  1  *  Fig.  12,  Attenuation  Coefficients  f 

wavelength  in  microns  Continental  Haze  Models  (From  Ref 
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begin  co  grow  rapidly  from  below  0.1  p  into  the  0.1  -  1  pt  range. 

There  is  no  one-to-one  relationship  between  relative  humidity  and 
meteorological  range,  but  an  approximate  upper  limit  on  aerosol 
attenuation  can  bo  deduced  from  measurements  by  Rosenburg  (Ref  30:5). 
From  his  measurements  of  horizontal  transparency  and  relative  humidity 
(Figure  13)  the  approximate  limit  on  minimum  true  visibility  due  to 
water  aerosols  can  be  determined,  as  in  Table  III. 
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vs.  Relative  Humidity  (Ref  30:5) 


Table  III 


Approximate  Limit  on  Visibility  vs.  Relative  Humidity 


%  Relative  1  emidity 

30 

40 

50 

60 

70 

80 

90 

100 

Maximum  cf  (Km  *) 

.05 

.15 

.3 

.8 

.8 

1.2 

1.8 

Minimum  o  (Km  ) 

.05 

.05 

.05 

.05 

.05 

.05 

.05 

,05 

Approx,  minimum 
true  visibility 
due  to  scattering 
(Km) 

80 

26 

13 

5 

5 

3.2 

2.1 

0 

It  is  in  the  region  of  80  -  100  %  relative  humidity  that  visibilities 


can  go  below  3  Km  and  the  particle  size  distribution  start  to  shift  its 
peak  from  below  0.1  p  to  above  2  p  for  clouds  and  fogs.  This  transition 
depends  on  the  ratio  of  soluble  to  insoluble  condensation  nuclei  and  tine 
recent  history  of  the  aerosol.  For  ire  not  within  the  t,tatc  of 
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the  art  at  this  time  or  in  the  forsceable  future. 

Therefore  a  linear  extrapolation  of  the  attenuation  coefficients  at 
0.42  ji  and  0.69  p,  vhen  corrected  for  Rayleigh  scattering  and  water 
vapor  absorption,  will  yield: 

1.  a  good  approximation  of  the  attenuation  coefficient  at  1.065  p,  for 
visibilities  in  haze  above  3  Km,  with  relative  humidities  below  90% j 

2.  a  good  approximation  in  fogs  or  clouds  or  haze  with  a  visibility 
under  1  Km  and  humidities  of  100%} 

3.  a  fair  approximation  in  the  transition  range  between  haze  and  fog, 
between  relative  humidities  of  90-100%,  and  between  visibilities  of 

1-3  Km.  In  this  region  the  extrapolation  will  slightly  overestimate  the 
attenuation  coefficient  at  1.065  p.  Duo  to  the  errors  inherent  in 
measuring  meteorological  range,  this  effect  is  probably  negligible. 

This  linear  extrapolation  method  yields  the  attenuation  coefficient 
at  one  ground  position.  This  can  be  inserted  into  the  bottom  of  the 
homogeneous,  exponential,  or  stratified  model  atmosphere  to  get  an 
attenuation  coefficient  as  a  function  of  altitude.  From  the  illumi.na.or~ 
target  range  and  target-receiver  range  this  can  be  used  to  yield  a  total 
atmospheric  attenuation.  This  will  be  accomplished  in  Chapter  4. 
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IV.  Lock-on  Range  in  a.  Real  Atmosphere 


Nominal  Lock-on  Ranges 


In  a  real  (attenuating)  atmosphere  the  transmission  through  the 
illuminator-target  optical  path  R  is  given  by: 


R 

T  =  10-*1/  S(R)  dR 

i  0 


(39) 


and  the  transmission  through  the  target-detector  path  R^  is  given  by: 
R1 

(AO) 


T„  .  dEd 

d  0 


The  clear~air  lock-on  range  R^c  given  in  Eqs  (20)  and  (23)  are  related 


to  the  true  lock-on  range  R  by: 

KdC  •  *dc\/v7 


(Al) 


The  system  parameter  Rdc  is  independent  of  weather.  The  weather 
parameter  Td  is  independent  of  receiver  system  design  because  the 
receiver  area  is  small.  is  very  slightly  dependent  on  target  area 
due  to  beamspreading,  significant  only  for  narrow  beams  illuminating 
large  targets. 

If  the  receiver  and  illuminator  are  together  in  the  same  aircraft. 


R  =  Rd,  Tl  =  Td  =  T,  so: 
Rdt  =  RdcT 


(A2) 


For  a  homogeneous  atmosphere,  often  encountered  below  a  strong 
thermal  inversion,  3  =  constant  and  Eq  (A2)  becomes: 

V  = 


(A3) 


for  all  dive  angles.  The  solution  to  this  equation  is  plotted  in  Figure 
1A  for  {3  from  0  to  AO  dB/Ku,  which  covers  Che  entire  range  of  flying 
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weather  conditions  normally  encountered.  At  short  ranges  there  will  be 
a  significant  but  not  prohibitive  differences  between  IT.  and  Rdc  for 
most  laser-guided  systems  for  most  weather  conditions.  At  ranges  beyond 
3  Km,  weather  conditions  can  easily  be  encountered  that  do  not  prohibit 
flying  by  visual  references  but  do  prohibit  a  reasonably  priced  laser- 
guided  system  from  locking  on  to  the  target  prior  to  launch. 

The  problem  can  be  worked  in  reverse,  going  from  a  statistical 
climatological  analysis  through  attenuation  coefficients  and  their 
statistical  variation,  through  this  chart  to  a  ’desired  clear-air  lock-on 
range  from  which  a  system  can  be  designed.  It  was  stated  once  before  in 
the  development  of  the  2ero-attenuation  case,  that  increasing  the  receiver 
sensitivity  would  increase  the  zero-attenuation  lock-on  range  by  the 
square  root  of  the  increase.  This  is  true  for  the  "clear  air"  case,  but 
for  the  attenuating  atmosphere  case  the  actual  increase  will  be  somewhat 
less.  A  knowledge  of  the  climatological  statistics  is  necessary  to 
predict  the  exact  amount.  If,  for  example,  from  Figure  14  the  expected 
atmospheric  attenuation  is  5  dB/Km  (about  2  mile  visibility)  and  a 
system  is  designed  for  a  zero-attenuation  lock-on  range  of  2  Km,  then  it 
will  have  an  actual  lock-on  range  of  0.8  Km.  If  the  receiver 
sensitivity  were  increased  by  100  times  (very  easy  to  do  for  present 
operating  systems)  then  the  zero-attenuation  lock-on  range  would  be 
increased  to  20  Km.  Following  the  same  5  dB/Km  attenuation  line  up  to  a 
new  zero-attenuation  lock-on  range  of  20  Km,  however,  yields  an  actual 
lock-on  range  of  2  Km,  or  what  the  system  was  designed  for  in  the 
absence  of  an  atmosphere.  If  the  visibility  wore  to  double,  the  actual 
lock-on  range  would  increase  by  65.^  for  this  combination  of  atmospheric 
attenuation,  system  lock-on  range,  and  homogeneous  atmosphere. 
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For  an  exponential,  unstratified  atmosphere,  P  =  P(h=0)eXp("h/1*2) 
(Ref  37s 7-1),  where  the  height  scale  factor  1.2  is  an  average  of  a  large 
number  of  field  measurements.  Therefore  if  0  is  the  aircraft  dive  angle 

and  the  height  h  =  R  sin  6,  then: 

P'dt  rd!Rsin_0 

R  =  R  ■10"'1f(R)dR  =  P.  lO"’1  P<h=0)J  e  dR  (44) 

dt  dc  0  dc  ^  N  ' 

Eq  (44)  was  solved  and  plotted  in  Figure  15  for  a  60°  dive  angle  and  for 
®(h=0)  ^rom  0  t0  ^0  dB/Km.  Plots  for  10°, 15°, 30°,  and  45°  dive  angles 
are  in  Appendix  C.  AFCRL  data  (Ref  37:7-3)  shows  that  this  model  can  be 
used  up  to  a  10  Km  altitude. 

The  graph  for  the  10°  dive  angle  is  nearly  the  same  as  the  graph 
for  the  homogeneous  atmosphere  case,  as  would  normally  be  expected;  the 
differences  that  occur  are  significant  only  at  the  longer  ranges  and 
middle  attenuation  coefficients.  The  graphs  shift  slowly  as  dive  angle 
is  increased  to  60°,  showing  steadily  improving  performance  as  the  dive 
angle  increases.  It  is  not  necessary  to  include  curvature  of  the  earth, 
as  some  have  suggested,  as  the  aircraft  is  not  physically  at  a  100  Km 
range.  If  the  system  is  designed  so  that  it  could  acquire  and  track  the 
target  at  a  100  Km  range  in  a  zero-attenuation  situation,  then  in  a  real 
atmosphere  with  the  aircraft  at  ranges  less  than  10  Km  the  system  would 
lock  onto  the  target  at  a  range  of  "x"  Km. 

The  improvement  of  lock-on  range  with  increased  dive  angle  is 
dependent  on  the  target  reflective  cross  section  not  increasing  as  dive 
angle  increases.  The  importance  of  this  fact  cannot  bo  overemphasized. 

It  would  be  very  misleading  to  imply  that  an  increased  dive  angle 
improves  lock-on  range,  even  if  all  atmospheric  conditions  are  constant. 
The  reflective  cross  section  of  nic  real  target  as  a  function  of  dive 
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angle  must  be  known  to  work  the  total  problem.  This  is  closely  related 
to  the  recommendation  on  page  19  for  follow-on  studies  of  real-target 
reflective  cross  section  analysis.  If  the  target  reflective  cross 
section  were  doubled,  then  the  problem  could  be  considered  as  if  the 
receiver  sensitivity  were  doubled,  which  was  discussed  on  page  49.  It 
is  recommended  that  the  results  of  the  follow-on  study  suggested  above 
be  inserted  into  this  sub-problem  and  that  charts  be  made  of  the  total 
problem  for  different  dive  angles  as  was  done  here,  but  including  the 
effects  of  the  changing  real-target  reflective'  cross  section  as  dive 
angle  is  increased. 

For  a  real,  stratified  atmosphere,  the  vortical  profile  of  a  test 
flight  over  Memmingen,  Germany  on  3  June  1970  (kef  18)  was  closely 
modeled  by} 

P  =  .01  dB/Km  above  5  Km  altitude 

=  .09  -dB/Km  between  1.7  -  5  Km  altitude  (45) 

(3  =  varied  in  stops  for  each  computation  but  constant  in 

each  computation  from  ground  level  up  to  1.7  Km. 

which  were  approximated  from  the  measured  total-scattering  functions 

for  red  (.68  p)  and  blue  (.48  ^i)  filters  on  the  test  aircraft.  On  this 

day  the  thermal  inversion  was  very  strong  and  located  at  an  altitude  of 

1.7  Km.  The  flight  was  in  the  afternoon  so  there  was  considerable 

thermal  convective  mixing,  as  evidenced  by  the  constant  scattering 

functions  below  the  inversion.  For  these  computations  the  bottom  layer 

attenuation  coefficient  was  set  at  various  values  to  create  a  family  of 

curves  for  this  typo  of  atmosphere.  If  the  aircraft  is  below  1.7  Km 

altitude,  a  homogeneous  atmosphere  results,  and  the  attenuation  is  =  |3R, 

for  h  =  R  sin  8  <  1.7  Km.  If  the  aircraft  is  between  1.7  -  5  Km  altitude, 

a  two-layor  stratified  atmosphere  results,  and j 
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hi  R 

sin  0  dt 

A  =  j31dR+Jg.  dR=  +  ,09Rdtsin  6  -  .153)  (46) 

0  1.7 

sin  3 

If  the  aircraft  is  above  5  Km  altitude,  a  three-layer  stratified 
atmosphere  results,  and: 

hi . 1-  P 

sin  0  sin  6  dt 

A  =  Jp1dR+Jg;)dR+  J*P3  dR  =  — (1. 731  +  .247  +  .OlR^sin  0)  (47) 

0  1.7  ~  5 

sin  0  sin  0 

This  stratified  atmosphere  problem  was  solved  and  plotted  for  a  60°  dive 
angle  in  Figure  16.  The  charts  for  10°,  15°,  30°,  45°  dive  angles  are 
in  Appendix  C.  The  second  inversion  at  5  Km  altitude  has  little  effect 
on  the  lock-on  range,  but  the  inversion  at  1.7  Km  is  very  conspicuous. 
Above  the  ma^n  inversion,  an  increase  in  system-desigr.ed  (zero- 
attenuation)  lock-on  range  has  an  almost  equal  increase  in  the  actual 
lock-on  range. 

This  form  of  computation  is  valid  for  any  number  of  inversions  at 
any  altitude.  The  three  cases  can  bo  combined s  an  exponential  atmosphere 
above  an  inversion  and  a  homogeneous  atmosphere  below.  Further,  if  it  is 
known  that  a  little  mixing  is  occuring  below  an  inversion  of  unknown 
strength,  the  homogeneous  atmosphere  chart  can  be  used  for  a  worst  case 
and  the  exponential  atmosphere  charts  can  be  used  for  a  best  case. 
Variation  in  Lock-on  Ranges 

The  nominal  lock-on  range  will  vary  considerably  J  roi.i  second  to 
second  in  the  tactical  combat  situation.  The  reflective  cross  section 
of  a  TSAC  jeep  model  with  an  olive  drab  paint  has  a  standard  deviation 
of  +  28% ,  so  the  lock-on  range  will  have  a  standard  deviation  of  14%  for 
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that  factor  alone,  even  neglecting  the  bright  specular  reflections.  In 
their  experimental  validation  of  their  model,  TSAC  (Ref  2s49,53) 
encountered  a  10%  fluctuation  in  the  power  output  of  their  military 
1.065  j.i  laser  illuminator.  This  author  using  a  shoulder-fired  laser 
illuminator  encountered  fluctuations  of  25%  many  times  in  every  5-second 
firing  test.  A  standard  deviation  of  these  laser  power  readings  is  not 
available,  but  will  be  approximated  at  +8%.  This  will  cause  a  standard 
deviation  in  lock-on  range  of  +4%  for  this  effect  alone. 

Variations  in  detector  threshold  are  extremely  temperature- 
dependent,  and  unless  temperature-compensated  electronics  are  employed, 
will  cause  considerable  variation  in  lock-on  range  for  seeker-limited 
systems.  Background-limited  systems  are  not  similarly  affected;  but 
they  have  other  limitations  which  depend  on  the  design  of  their  circuit 
and  imaging  systems.  The  variation  in  cos  ©d  is  assumed  to  be  small, 
hopefully  due  to  proper  system  design.  If  the  angle  @d  were  off  by  an 
arbitrary  20°,  then  cos  6d  =  . 94,  resulting  in  a  maximum  of  3%  reduction 
in  zero-attenuation  lock-on  range. 

Atmospheric  fluctuations  fairly  well  average  out  when  the  atmos¬ 
pheric  visibility  is  over  3  Km,  but  may  fluctuate  locally  up  to  a  factor 
of  six  when  the  atmospheric  visibility  is  under  0.5  Km.  This  is  an 
unpredictable  but  very  significant  factor  when  attacking  a  ground  target 
surrounded  by  a  shallow  but  dense  ground  fog  or  a  thin  stratus  layer. 

The  biggest  variations  in  determining  the  attenuation  coefficients  ares 

1.  Estimating  the  meteorological  range  through  optical  filters.  This  can 
easily  cause  an  under-  or  overestimation  in  the  attenuation  coefficient 
by  a  factor  of  two  at  the  observer’s  position. 

2.  Identifying  the  vertical  structure  of  the  aerosol  concentration  of 
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the  atmosphere,  from  the  observer  up  to  at  least  above  the  first  strong 
inversion  layer  or  stratus  layer. 

This  uncertainty  in  the  weather  element  car  be  greatly  reduced  by 
bouncing  a  lidar  radar  at  1.065  p  off  a  large  tethered  weather  balloon 
as  it  is  raised  through  the  bottom  layer  of  the  atmosphere,  and  measuring 
the  real  attenuation  in  that  atmosphere  as  a  function  of  altitude.  The 
reflectivity  of  the  balloon  could  be  increased  to  nearly  80%  by  spraying 
with  aluminum  paint.  As  the  balloon  is  spherical,  stabilization  is 
unnecessary.  Since  the  balloon  is  tethered,  it  can  be  recovered.  This 
method  should  be  very  successful  when  the  first  strong  inversion  is 
below  1000  feet  (0.3  Km),  as  it  often  is  in  central  Europe.  It  was 
intended  to  experimentally  test  this  concept  but  equipment  and  weather 
problems  prevented  actual  tests.  The  experiments  actually  conducted 
are  described  in  the  next  chapter. 
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V.  Experimental  Test:  of  Real  Atmosphere  Problem 

Experimental  Concepts  to  be  Tes  ted 

Goal  it  I  was  tc  measure  the  ground-level,  horizontal  atmospheric 
attenuation  at  1.065  ^i,  0.695  fi,  and  0.42  ;j.  simultaneously  in  different 
Jiaze  and  fog  conditions,  to  determine  if  attenuation  at  1.065  ft  could 
be  determined  from  measurements  of  meteorological  range  through  narrow- 
band  filters  in  the  visual  portion  of  the  spectrum.  Attenuation  at  the 
first  wavelength  was  to  be  measured  directly  by  firing  a  shoulder- 
mounted  laser  illuminator  at  a  receiver  from  three  known  distances. 

Goal  it2  was  to  measure  the  slant-range  atmospheric  attenuation  at 
1.065  fi  by  bouncing  the  illuminator  beam  off  an  aluminized  weather 
balloon  and  analyzing  the  returned  signal  with  a  tolephotometer  and 
oscilloscope.  Due  to  the  lack  of  funds  and  time,  suitable  receiver 
equipment  could  not  be  obtained  during  the  short  thesis  period,  so  goal 
it2  could  not  be  achieved. 

Experimental  Organization  and  Procedure 

The  shoulder-fired,  NdsYAG,  repetitively  q-switched  laser 
illuminator  had  a  diffuso"  cap  which  spread  the  laser  beam  over  a  60° 
cone.  This  illuminator  was  fired  at  the  receiver  from  distances  of  70, 
200,  and  300  feet.  The  power  was  received  by  an  EGSC  spectroradiom^ter 
with  a  100  j?  monachromator  centered  on  the  laser  line  at  1.065  +  0.002|i. 
The  receiver  had  a  14°  field  of  view  (0.18  steradian  solid  angle)  and 
was  eloctrocryogenicaily  cooled,  with  an  S-l  spectral  response.  Received 
signal  strength  was  read  or  the  large  internal  meter;  it  was  not  possible 
to  modify  the  borrowed  equipment  for  any  other  output. 

Measurements  were  made  ori  the  roof  of  the  AFIT  School  of  Engineer ii.t* 
Building  640,  Wright-Patterson  AFB,  elevation  1004  feet  above  rxan  sea 
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level.  Only  three  days  between  1  Oct  1971  and  7  March  1972  had  suitable 
weather  for  attenuation  measurements  in  haze:  26  Oct  1971,  10  Jan  1972, 
and  20  Jan  1972.  borrowed  equipment  was  not  available  on  the  first  two 
days  and  from  2/i  Jan  -  7  March  1972.  Therefore  only  one  day  was 
available  for  measurements,  although  equipment  ard  procedures  were 
tested  on  several  ether  days.  On  20  January  the  white  light  visibility 
varied  from  800  feet  in  dense,  warm,  postfrontal  fog,  to  7  miles  in 
light  haze  below  a  thick  overcast  stratus  layer.  The  temperature  varied 
from  51°“5A°F  and  the  relative  humidity  varied ‘from  93%-100%. 

A  spectral  scan  of  the  background  from  1. 0-1.1  p.  was  made  at  1200 
hours  on  a  heading  of  085°  (magnetic)  and  at  an  uplook  angle  of  5°  to 
clear  the  horizon,  into  a  moderate  fog  fusing  into  a  very  dark,  thick, 
low  overcast. 

Five  complete  sots  of  measurements  were  made.  In  each  set: 

1.  The  meteorological  ranges  in  white  light,  through  narrowtand  red, 
narrowband  blue,  and  wideband  red  filters  were  measured  directly  by  the 
author.  The  distance  to  the  farthest  visible  objects  was  determined 
from  base  engineering  maps  and  local  aeronautical  charts.  Several 
distant  buildings,  telephone  poles,  tree  lines,  or  other  black/white 
objects  were  used  to  determine  observed  (not  the  true)  meteorological 
range.  If  time  and  budget  had  permitted,  a  sot  of  large  black  and  .date 
panels  would  have  been  sot  up  in  a  line  and  photographed  during  the 
experiment  through  the  rod  and  blue  filters  for  a  later  microdensitometer 
analysis  of  contrast  degradation.  The  intent  of  the  experiment  was, 
however,  to  determine  whether  estimates  of  the  visibility  through  rod 
and  blue  filters  could  give  an  approximate  attenuation  coefficient,  by 
quick  ar.d  simple  field  meteorological  measurements. 
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2o  The  laser  illuminator  with  the  diffusor  cap  on  was  fired  at  the 
receiver  from  distances  of  70,  200,  and  300  feet;  the  received  signal 
strength  was  recorded.  Background  strength  was  recorded  before  and 
after  each  laser  firing.  Unfortunately  there  was  a  25%  random 
fluctuation  in  the  laser  signal  (This  had  not  occurred  during  prior 
trial  tests.),  and  the  meter  needle  would  fluctuate  randomly  within  a 
25%  range  many  times  per  second.  It  was  necessary  to  "eyeball  average" 
over  a  5-second  interval.  Either  the  laser  was  not  firing  every  time, 
or  the  pulses  were  not  uniform  in  power  and  total  energy.  The  meter 
maximum,  minimum,  average,  and  uncertainty  of  the  average  were  recorded. 
3.  A  detailed  description  of  the  weather,  temperature,  and  relative 
humidity  were  recorded.  These  were  compared  with  the  Wright-Patterson 
official  weather  observations  taken  5  Km  northeast,  and  with  the 
vertical  temperature  profiles  taken  17  Km  Northwest  at  Dayton  Municipal 
Airport. 

The  receiver  was  scanned  +5°  across  the  transmitter  with  no  change 
in  readings.  The  transmitter  was  scanned  +5°  across  the  receiver  with 
no  change  in  readings.  Therefore  aiming  of  the  transmitter  or  receiver 
was  not  a  problem. 

Analysis  of  Data 

The  spectral  scan  of  the  background  sky  is  presented  in  Figure  17, 
which  appears  to  be  a  combination  of  the  severely  attenuated  solar 
influx  which  decreases  with  wavelength  in  this  region,  a:id  the  black- 
body  radiation  of  the  overcast  cloud  layers  which  increases  with 
wavelength  in  this  region. 
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Fig.  17.  Spectral  Scan  of  the  Background  Sky 


The  experimental  data  is  tabulated  in  Table  IV.  The  attenuation 
coefficients  are  determined  from  Eq  (38),  for  the  wavelengths  in  the 
optical  region.  The  laser  beam  attenuation  coefficients  can  be 
determined  by  (from  Appendix  D): 


in  dB/units  of  R  (48) 


where: 

1^  =  intensity  (irradiar.ee)  at  range  R^ 
-  intensity  (irradiance')  at  range  R^ 
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was  not  available  and  the  weather  was  not  sufficiently  hazy.  Since  only 
the  relative  intensities  for  the  three  distances  is  used  in  calculations, 
and  if  it  can  be  assumed  that  the  problem  with  the  laser  illuminator  was 
constant  for  all  the  firings  (as  it  seemed  to  be)  then  the  average  of 
the  meter  fluctuations  is  meaningful. 

The  vortical  temperature  profile  of  the  atmosphere, taken  at  0700 
and  1900  EST,  20  Jan  1972,  and  0700  EST,  21  Jan  1972,  from  Dayton 
Municipal  Airport,  appears  in  Figure  19.  A  warm  front  passed  Dayton 
about  1000  EST  on  20  Jan  1972. 

The  spectral  response  of  the  three  filters  used  in  this  experiment 
was  measured  on  a  spectral  transmissometer,  pointwise  corrected  for  the 
zero  and  100%  lines,  and  reconstructed  in  Figure  20.  The  narrowband 
blue  filter  greatly  reduced  the  apparent  brightness  of  the  surrounding 
landscape  and  made  visibility  judgement  difficult.  The  spectral  response 
of  the  blue  and  violet  landscape  also  made  visibility  judgement  difficult. 
The  narrowband  red  filter  likewise  had  the  same  problem  as  the  blue  and 
was  too  far  towards  the  infrared.  The  wideband  red  filter  had  excellent 
characteristics;  determination  of  meteorological  range  was  far  more 
precise  and  much  easier.  Measurement  of  visibility  through  wideband  red 
(0.60-0.70  p)  and  wideband  blue  (0.40-0.50  p)  filters  will  probably  yield 
much  more  accurate  results  than  the  narrowband  filters  selected  for  this 
experiment.  Adverse  effects  of  the  spectral  characteristics  of  the 
surrounding  visual  references  should  also  be  reduced  with  the  wideband 
filters. 

Results  of  Test 

Due  to  weather  and  equipment  problems  the  results  of  this  test  are 
inconclusive.  The  meteorological  range  through  the  blue  filter  was  too 
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Elevation  of  test  site  1004  MSI,* 

Elevation  of  Dayton  Muni*  1003  MSL. 
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Temperature  in  Degrees  C. 


Fig.  19.  Atmospheric  Vertical  Profile 
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low  to  match  the  theory.  The  spectral  response  and  inherent  com:  a. 

through  red  and  blue  filters  of  the  visual  references  around  the  te^; 

site  were  not  measured  due  to  a  lack  of  time.  Additionally,  it  was 

desired  to  develop  a  system  applicable  to  general  field  use  by  uncra:  - 

personnel  and  with  an  absolute  minumum  of  equipment. 

This  test  needs  to  be  repeated  with  a  stable,  single-mode,  l.0. 3 

laser  illuminator,  and  wideband  red  and  blue  optical  filters.  A  CW  1 

is  recommended.  /  Reproduced  frorrT  IsS 

'  LLcs'  *vailab!nJ:rY  § 


Fig.  20.  Pointwise  Reconstructed  Spectral  Scan  of  Filters  L'^ed 
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VI.  Conclusions  and  Recommendations 

1.  The  most  serious  reduction  of  the  lock-on  range  at  1.065  p.  is  caused 
by  the  aerosol  attenuation  of  haze,  fogs,  and  clouds.  The  most 
important  elements  of  aerosol  attenuation  are  the  aerosol  composition 
and  size  distribution,  which  vary  extremely  with  time,  location,  and 
altitude.  Determination  of  aerosol  attenuation  coefficients  from 
routine  nonvisual  meteorological  measurements  is  beyond  the  state  of  the 
art.  Relative  humidity  is  not  a  useful  indicator  of  aerosol  composition 
or  attenuation.  It  should  be  possible  to  determine  attenuation 
coefficients  from  an  extrapolation  of  the  coefficients  at  nearby 
wavelengths  in  the  visual  range,  which  could  be  determined  from  measuring 
the  visibility  through  wideband  red  and  blue  filters.  This  concept  was 
tested  against  calculations  and  measurements  of  other  researchers  with 
successful  results  for  most  cases.  Experimental  tests  of  this  concept 
were  inconclusive  due  to  equipment  problems  and  a  lack  of  suitable 
weather.  It  is  recommended  that  the  experiment  be  repeated  with  wideband 
filters  and  a  stable,  CW,  laser  illuminator- 

2.  Lock-on  ranges  are  strongly  affected  by  changes  in  the  target 
reflective  cross  section,  which  may  change  by  a  factor  of  2  -  3  times 
for  different  illumination  and  observation  angles.  The  reflective  cross 
sections  of  dirty  or  v  „  surfaces  are  unpredictable  at  present.  It  is 
recommended  that  the  reflective  cross  sections  be  computed  for  several 
enemy  threat  vehicles  for  all  viewing  angles,  for  clean,  wet,  and  dirty 
surfaces.  It  is  recommended  that  the  reflective  cross  section  as  a 
function  of  dive  angle  be  determined,  and  this  function  be  inserted  in 
the  total  atmosphere  lcck-on  range  problem  of  Chapter  IV. 

3.  The  vertical  structure  of  a  i\ul  atmosphere  can  bo  bracketed  between 
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an  exponential,  unstratified  atmosphere  (for  no  thermal  mixing  or 
convective  activity  for  several  days)  and  a  homogeneous  atmosphere 
(considerable  convective  activity  below  a  strong  thermal  inversion). 
Computer  solutions  for  both  ..ases,  and  for  a  measured  atmosphere  from 
flight  test  data,  were  computed  for  five  dive  angles  between  10°-  60° 
and  for  surface-measured  attenuation  coefficients  between  0—^0  dB/Km. 
Inhomogeneities  in  the  atmosphere  average  out  when  the  visibility  is 
over  3  Km,  but  can  cause  fluctuations  of  25%  below  2  Km  visibility  and 
fluctuations  of  600%  below  0.5  Km  visibility.  • 

4.  The  lock-on  range  can  be  increased  several  orders  of  magnitude  by 
using  cryogenically  cooled  optical  arrays  and  electronic  circuits  that 
discriminate  against  extended  area  returns  and  false  targets  from  the 
highly  reflective  background  or  foreground  objects  also  illuminated  by 
the  laser  beam. 

5.  The  contrast  threshold  characteristics  of  the  photopic  eye  as  a 
function  of  wavelength  are  not  known.  It  is  recommended  that  a  black/ 
white  (black  velvet/white  magnesium  oxide)  target  be  viewed  through  many 
narrowband  and  wideband  filters  to  determine  the  contrast  threshold 
characteristics  of  the  photopic  eye  as  a  function  of  wavelength  and 


filter  bandwidth. 
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Appendix  A_ 

Comparison  of  Beam  Divergence  Definitions 

A  computer  program  was  developed  to  integrate  the  energy  of  a 
gaussian  beam  th?c  '..'as  contained  within  a  cone  of  increasing  angular 
dimension  centered  in  that  beam.  If  K  is  the  portion  of  the  beam’s 
energy  within  the  cone,  then  (from  Eq  (2)): 


K  = 


2it  r* 
c_ 

-  w  frww 
0  0 


Q~c('u'>  rdrdG  =  f  ~  e~cy  dr  =  1  -  e"C(w^ 

J  w“ 

0 


(50) 


At  the  90Z  one losed-energy  points,  K  =  .9000,  so  when  -  =  \  c  =  In  10 

w 

=  2.3026.  The  integrating  program  was  developed  at  this  point  to  test 
the  program  against  exact  integral  results  just  performed.  The  program 
integrated  Eq  (50.)  by  a  gauss  15-point  integration  subroutine  GLDI  from 
the  AFIT  computer  center,  with  an  interval  size  1/5  of  the  integration 
limits,  thus  creating  an  effective  gauss  75-point  integration.  The 
values  if  K  were  computed  from  0  <  £  <  2.5  in  steps  of  0.05  by 
repetitive  iteration.  These  values  agreed  with  the  exact  solution  of 
Eq  (50)  to  the  third  decimal  place.  The  program  was  used  soon  after¬ 
wards  for  difficult  integrals  that  could  not  bo  checked  against  the 
exact  solutions:  a  gaussian  beam  illuminating  a  rectangular  target. 

Now  to  compare  the  relative  widths  of  the  gaussian  curves.  For  a 
beam  defined  by  Eq  (2),  with  r  as  the  only  variable,  losing  at  the 
gaussian  exponential  ar.d  therefore  neglecting  the  integrating  height- 
normalizing  constant  when  comparing  beam  widths: 

i Or)  _  — 2.3026(~)  vj  j  w  •« 

Ho)  -  1  "  ><"'?=  P..3020/  !„  (51) 

Eq  (j1)  gives  the  relationship  between  column  2  and  coluui  5  in  Table 
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9 

I.  For  example,  row  5:  by  this  definition  I(r)/I(o)  =  l/e“  =  0.1353. 
Therefore: 

-=  2.3026/Vln  (e2)  =1.075 

r 

which  is  the  multiplicative  factor  for  converting  this  defined  laser 
beam  width  to  the  standard  definition  of  beam  width  for  this  thesis. 
By  setting  r  =  w  and  letting  I(r)/I(o)  be  as  defined  for  each  row. 


I(r)  _  -c 
Ko)  '  e 


,  or 


c  =  In 


I(o) 

Kr) 


(52) 


Eq  (52)  gives  the  relation  between  column  2  and  column  3  of  Table  I. 

2 

For  example,  row  5.  By  this  definition:  ,t(r)/I(o)  =  1/e  .  Therefore 
c  =  2.0000.  Finally,  for  each  row,  setting  r  =  w,  and  having  computed 
c  from  Eq  (52),  Eq  (50)  can  be  used  to  determine  the  total  energy 
enclosed  by  the  cone,  in  column  4  of  Table  I. 
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Appendix  ji 

Integration  of  Beans  Across  Targets  of  Various  Shapes 

For  a  constant-irradiance  cone  of  light  falling  on  the  center  of  a 
circular  target,  the  ratio  of  power  intercepted  by  the  target  (P  )  to 

transmitted  power  (P)  (neglecting  any  atmospheric  effects)  is  given  by: 
•  .  2 

K  =  p—  =  Or)  for  h  <  d 

(53) 

K  s  1.  for  h  >  d 

where  h  =  the  vertical  dimension  (height)  of  the  round  target,  and  d  = 
the  beam  diameter  at  the  target  position. 

For  a  constant-irradiance  beam  of  diameter  d  illuminating  a 
rectangular  target  of  height  h  and  length  nh.  There  are  four  possible 
cases,  diagramed  in  Figure  21. 


& 
2 


Fig.  21.  Geometry  of  Constant  Beam  on  Rectangular  Target 
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Case  a.  When  h  >  cl,  all  the  beam  energy  is  on  the  target,  and  K  =  1. 

T.  li 

Case  b.  When  ^<~<  1,  then  for  each  quadrant,  all  the  beam  energy 
between  0  <  6  <  6,  is  intercepted  plus  the  energy  in  the  triangle 
between  <  9  <  90°.  If  the  total  beam  energy  is  normalized  to  unity, 
then  the  irradiance  is  H  =  4/ud*"#  sinc0  ^  =  arcsin(-). 


-J H  dS  =  Al  t~~2  rclrc^  +  l~ 

S  J  J'l Td“  J  h  Jn 

A  A  A _ 


r=h/2  sin  9 


r  dr  d0 


02=arcsin(— )  0 


2  .  .h.  .  2h 

—  arcsin(-rj  +  —r 
tt  d  -rrd 


A  -  <§>' 


Case  c.  When 


X  1\  X 

•====  <  -7  <  —  »  then  the  energy  intercepted  in  each 
1  +  n2 


quadrant  is  that  in  the  triangle  from  0  <  6  <  6  ,  plus  the  circular 
sector  from  6,  <  6  <  02,  plus  the  triangle  from  @2  <  0  <  90°.  Since 
0X  =  arccos(^p). 


01=arccos(“)  r=nh/2  cos  6  e2=arcsiri(^)  r=cl/2 


K  =  4 


Py* 

JiTd 


d9  +  4 


S.^arccosC—) 
1  d 


•dr  d0 


r=h/sin  9 


— 9  dr  d0 


02=arcsiu(^-) 


„  2  ' 
7T 


tircsin(^)  -  arccos(— ) 


nh]  +  2nh 

J  7TC 


+Sr\/i-(T) 


Case  d.  When  all  the  target  is  in  the  Learn,  the  ratio  of  target  area 
2  9 

(nh  )  to  the  beam  area  (-rrd“/4)  is  the  fraction  of  energy  intercepted, 

K  -  f  (|)2  (50) 

For  a  gaussian  beam  on  targets  of  different  shapes,  the  total  power 


intercepted  by  the  target  can  b'  solved  by  one  cafe,  since  the  gain  c’ an 
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beam  has  no  discontinuities  and  is  very  much  larger  than  the  target 
so  can  be  considered  filling  the  transmitted  hemisphere.  For  the  round 
targets  the  problem  is  the  same  as  in  Appendix  A: 

K  -  1  -  (57: 

For  the  rectangular  targets,  no  exact  solution  exists,  so  a  numerical 
integration  (from  Appendix  A)  was  used  at  each  of  100  points  from  n  = 

0  to  n  =  2.5,  as  for  the  constant-irradiance  problem  above.  Using  the 
same  nomenclature  as  above,  with  0^  =  at  the  corner  of  the 


rectangle,  so  0^  =  arccot  n, 

0  =arccot(n)  r=hn/cos(0) 
i  /->  n 


K  =  4 


0  =arcuot(n) 

1  r 


[i'll 

J  lid2 


n/2  r=h/sin(0) 


2104r  -9. 2104  (~)  .  ^  , 

—7 - e  d  drd0+  4 


if 

J  .  Jinc 


integrand 


dr  d0 


01=arccot(n)  0 


=  1  - 


2.3026 


"9.21  Ifni  17 yr  p 3. 2104  h  j 

~  d2  cos“0. _  1 2. 302o_  e~  ?  siirej, 


0  0^=arccot(n) 

A  computer  program  was  created  to  solve  and  plot  Eqs  (54)  -  (58)  for 
values  of  n  (where  appropriate)  from  0  to  2.5  in  steps  of  0.025,  for 


a  circle,  square,  and  rectangles. 
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These  eight  charts  are  discussed  on  pages  51-53. 
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Computation  of  Attenuation  Coefficients 


If  1^  =  intensity  received  at  near  distance  R^,  I2  =  intensity 
received  at  far  distance  R7,  and  3  =  attenuation  coefficient  in  dB/unic 
length  of  R^  &  R2  in  the  experiment  in  Chapter  V,  then: 

1  riy  ,  jpy 

I,  =  10  L  10J  ,  or  I  =  IrJ  10  L  10J  (59) 


1  Ri 


I  _f^R2* 
=  ~  m~L  in- 


ft]. 


I-  =  -r  10"L  10-1  ,  or  I  =  I9R9  10+L  10 

2  2  0  2  2 

R2 


Therefore:  r  (pR  1 

I  R?  loHaJ  =  I^RplO  10  ,  or  10^(R2  "  Rl)/10  = 


10  m  ri 

^  (R2  —  Rx)  Iog10  I2  ii,  (*s; 

For  small  A3,  A^,  AI9:  A3  =  §f  A^  -  §§^Al2  (62] 

The  negative  sign  of  the  second  term  is  necessary  because  Aj^  and  Al2 
have  opposite  effects}  an  increase  in  I1  alone  will  increase  3,  but  an 
increase  in  I9  alone  will  decrease  3. 

The  relationship  between  3  and  visibility  is  as  follows.  The  true 
meteorological  range  (true  visibility)  is  the  range  at  which  contrast 
is  reduced  to  2%  (note  some  countries  use  5 %),  or: 


0.02  =  10“^'/ 10 


,  or  3V  =  10  log..  50  =  17.0  •  or  p  -  — 

1U  V 
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